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METHOD OP SPLIT RIGIDITIES AND ITS APPLICATION TO 
VARIOUS BUCKLING PROBLEMS 
By P. P. Bijlaard 

SUMMARY 


A comprehensive treatise on the method of split rigidities is pre- 
sented. First the principles upon which the method is based are discussed. 
It is shown on new examples how these principles are applied. These 
applications are divided into problems where all component modes into 
which the actual behavior of a composite structure is split have the 
same boundary conditions and into those where these boundary conditions 
differ. Examples of the first type include sandwich columns with various 
boundary conditions, columns with batten plates, and latticed columns j 
examples of elastic and plastic buckling of sandwich plates with ortho- 
tropic core and of corrugated- core sandwich plates are also given. This 
type includes problems based on the same principles where only one mode 
has to be considered. As an example, the buckling stress of homogeneous 
plates under nonhomogeneous stresses in their plane is expressed in terms 
of their critical stress under homogeneous compression. To this group 
also belongs the determination of the ultimate load of plates under 
compression. An explicit formula is derived for the buckling stress of 
stringer panels which is a new example of the second type of problem. 

The problems were chosen so that the correctness of the method, which 
is basically an approximate one, can be shown by comparison with exact 
calculations or tests. 


INTRODUCTION 


In several papers a method has been used for calculating the buckling 
stresses of structures that buckle in composite modes which is called the 
"method of split rigidities." The method consists of splitting the 
buckling deflection into two or more component modes and expressing the 
buckling stress in terms of the critical loads for these component modes. 
References 1 to 19 are based partly or completely on this method. 

For example, in a sandwich plate (refs. 9 > 11> 12, and l4) there 
are three different buckling modes: (l) That of the single faces. 
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(2) that of the complete sandwich plate, assuming the faces to have zero 
flexural rigidity but assuming the core to have infinite shear rigidity, 
and (3) that of the complete sandwich plate, assuming the faces to have 
zero flexural rigidity and assuming the faces to be infinitely rigid 
against extension. In the case of a T-section (ref. 13) the three modes 
are: (l) Bending about the Y-axis situated in the plane of the web, 

(2) twisting about the shear center axis, and (3) plate buckling. In 
the case of buckling by general instability of a long cylindrical shell 
with stiff ening rings under external pressure (ref. 19) the component 
modes are bu ckling as an orthotropic shell, assuming the rigidity of the 
rings to be uniformly distributed along the length of the shell, and 
bu ckling of the shell between the rings. In the case of stringer panels 
(ref. 18) practically exact solutions are obtained Icy this method for 
conventional panel buckling as well as for forced crippling (refs. 16, 

18, and 20). 

Formulae for buckling loads or stresses have been derived in earlier 
papers and new ones will be presented herein by establishing equations 
between external and internal actions. In references 1, 2, 3, and 6 the 
external arid internal actions considered were the external and internal 
bending moments. In references 7 and 13 bending as well as torsional 
moments were considered. Bending moments in a cross section as well as 
deflecting and restraining transverse forces acting on an element were 
considered in references 4, 9 , H, 13, and 16. The actions and reactions 
considered in references 14, 15, 17, and 19 were the amounts of work done 
by deflecting and restraining forces. Columns and plates with initial 
deflections were considered in a footnote of reference 1 and in 
reference 15. 

The method of split rigidities was also applied in the postbuckling 
range (refs. 15 and 17). It is very powerful and leads to simple results 
in cases where a solution by existing methods is practically impossible. 

It was successfully used in two other projects (refs. 13 and 17). Until 
the present only a few other investigators have used this method (refs. 21, 
22, and 23) probably because no general comprehensive treatise on it has 
yet been made available.- Therefore the purpose of the present report is 
to give such a general explanation, with examples of the 'method of appli- 
cation in various special cases. 

The present investigation was carried out at Cornell University 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. The examples dealing with the 
ultimate load of plates under compression and with buckling of stringer 
panels are part of earlier work done for the Bell Aircraft Corporation. 

The author wishes to express his appreciation for the valuable coopera- 
tion of all concerned and to the Bell Aircraft Corporation for its per- 
mission to publish the abovementioned examples. 
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total cross section of column 

plastic redaction coefficients for plastic buckling of plates 
cross section of diagonal 

cross section of stiffener 
cross section of vertical 

half wave length of buckling for plate 

width of plate; also, width of bay of stringer panel 

effective width of plate 

constant; also spring constant (foundation modulus) 

center- to- center spacing of batten plates in built-up column; 
also, free length of single struts in latticed column 

effective free length of single struts between batten plates 

deflecting force 

ratios of shear to shear angle 

elastic modulus 

elastic modulus of corrugation 

elastic modulus of faces (skin) in sandwich plate; also, 
secant modulus in section "Ultimate Load of Plates Under 
Compression" 

tangent modulus 

tangent modulus of corrugation 
modulus of rigidity of sandwich core 
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G x ,Gy moduli of rigidity in X- and Y-directions of sandwich core 

H total thickness of sandwich column or sandwich plate, t + 2h 

h center- to- center distance of single struts of built-up column; 

also, thickness, of single face of sandwich column or sandwich 
plate 

I moment of inertia 

I c moment of inertia of single strut, of built-up column 

If. moment of inertia of single face of sandwich plate 

*r moment of inertia of reduced built-up column, Ah 2 /4 


I 6 moment of inertia of reduced sandwich column or sandwich 

plate, (l/2)h(t + h)^ 

K constant 

k buckling stress coefficient defined by equation (l^2) 

k^k-j^kg buckling stress coefficients referring to cases 0, 1, and 2, 

respectively 

k 8 buckling stress coefficient defined by equation ( 1^5) 

k-i^jkp^, . . . coefficients in equation ( 104) 


L 

l 

l r /r 

l eq L / r 

M 

M x ,My 


effective length of column 

length of column; also, length of stringer panel 

effective slenderness ratio for reduced case of built-up 
column 

effective slenderness ratio of built-up column 
moment 

bending moments 


M X y,My X 


torsional moments 
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m * (W r )/(°/ rc ) 

N flexural rigidity of homogeneous plate; also, flexural rigidity 

of face of sandwich plate 

N s flexural rigidity of reduced sandwich plate, given "by 

equation (59) 

n = lj c 

n' number of hays of stringer panel 

P buckling load 

p,q coefficients in equations ( 108 ) and (203) 

Q shear force; also, transverse shear force in plate 

Q' transverse shear force that would cause a unit angular 

distortion in latticed column 

R restraining force 

r radius of gyration of reduced built-up column, h/2; 

r = r x = r y if r x = r y 

r c radius of gyration of single strut of built-up column 

r x ,r y ratios defined by equations ( 78 ) and (8l) 

? x ,f y ratios defined by equations (120) and ( 123 ) 

S spring 

S,S' E ^,S' EV quarter spring constants for plates 

t core thickness of sandwich column or sandwich plate; also, 

thickness of homogeneous plate 

Uj^ normal, function 

v translation of single struts per -unit shear force in 

batten plates 


w 


deflection 
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x,y,z coordinates 

X,Y directions 

Y function of y only 

a angle 

a^ctg coefficients defined "by equation (a6) 

a m coefficient in equation (102) 

3 factor in equation (35) 

T aspect ratio of buckle in buckled plate, a/b 

y angular distortion; also, ratio defined by equation (221) 

5 elongation 

e strain 

e m membrane strain 

£ value «1 in equation ( 228 ) 

?j reduction factor for plasticity for concentric buckling of 

plates 

r) reduction factor for plasticity for finite deflection of plates 

0 angle; 9 = 9 X = 0y if 9 X = 0y 

0 x'®y ratios defined by equation (75) 

k ratio defined by equation ( 226 ) for stringer panel 

A = m/a 

p ratio defined by equation (95) 

v Poisson's ratio 

\ slope of column in case 2 of latticed column 

a normal stress 
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*e 

normal stress at edge 

- 


membrane normal, stress 


a z 

normal stress in plane parallel to middle plane of sandwich 
plate 


T 

coefficient in equations ( 151 ) and ( 15 k) 


$(x>y) 

function of x and y 


<P 

ratio defined by equation (211) 


t 

defined by equation (212) 


n 

ratio defined by equation ( 225 ) for stringer panels 


Subscripts : 


a 

actual (fig. 10(b)) 


b 

bending 


cr 

critical 

* 

e 

external 


eq. 

equivalent 

M 

h 

homogeneous 


i 

internal 


m 

maximum 


mid 

center (middle) 


P 

postbuckling; also, refers to fictitious end deflections as 
shown in figure 10(b) 


r 

reduced case 


s 

skin (face) for sdndwich plate; also, secant in section 
"Ultimate Load of Plates Under Compression" 


st 

stiffener 


u 

uniform 


ult 

ultimate 

r • • 
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x,y refer to X- and Y-directions, respectively; also, refer to X- 

and Y- strips , respectively, for orthotropic sandwich plate 

ys yield stress 

0,l,2 ,...n refer to cases 0, 1, 2, . . . n, respectively 

DESCRIPTION OF METHOD 


The method of split rigidities as applied to "buckling problems may 
be described as follows: 

(1) Splitting into two or more component cases: The elastic or 

elastoplastic behavior of a composite structure is split int o two or 
more component cases for which the individual buckling stresses can be 
easily determined. For example, the deflection w of a sandwich plate 
with membrane faces is split into its deflections w^ from bending alone 

(case l) and w 2 from shear deformation alone (case 2). It is assumed 
that in the composite structure the component buckling deflections w-^ 
and w 2 have the same shape as in cases 1 and 2 where they occur alone. 
The accuracy of the method will be greater the better this assumption 
is fulfilled. If it is exactly fulfilled, which as explained in this 
report is only possible if w-^ and w g have the same shape, the method 
is exact. 

(2) External actions and internal reactions: With the shapes of 

the component deflections for the combined case thus determined, it is 
always possible to establish equations between the external actions and 
the internal reactions. These actions and reactions may be the external 
and internal moments acting in a cross section, the deflecting and 
restraining forces acting upon a small element, or the work done by 
deflecting and restraining forces. 

(3) External action* If the actual combined case is split into two 
or more, say n, component cases with deflections w-^, Wg, w^, . . . w n , 

the external action is proportional to the actual buckling load F cr for 
the combined case and may be expressed in terms of P cr and the indi- 
vidual deflections w-^, Wg, Vy . . . w q . For example, splitting the. 

deflection w of a sandwich column with membrane faces into the deflec- 
tion Wj_ from bending and the deflection Wg. from shear deformation, 
the external moment M g for the combined case is equal to P cr (Vq + Wg\ . 
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(4) Internal action: As stated under item (l), in the combined 

case with a deflection w = w^_ + Wg, the component deflections and 

Wg axe assumed to have the same shapes as in the separate cases 1 and 2 
with buckling loads P^ and Pg. The internal reactions depend only on 
the shape of the deflection and not on the magnitude of the compressive 
force. Therefore, the internal reactions in the combined case axe equal 
to those in case 1 alone since additional deformations (case 2) increase 
the deflections, but not the internal reactions. In the same way it 
follows that the internal reactions are equal to those from case 2. For 
example, in a sandwich column with membrane faces, the deflection Wg 

from shear deformation increases the deflection beyond w 1 from bending, 
but it does not increase the internal moment from w-^, and vice .versa, 
as shown extensively in reference 11. 

As will be explained later, the rigidities of the structure against 
the partial deflections w-^ and Wg can be considered to be supplied 

by sets of springs S-^ and Sg with different spring constants 
and Cg acting in series to restrain the deflection of a column without 
any proper flexural rigidity (fig. l). Hence it is obvious that the 
internal reaction for the combined case, that is, of the combined springs, 
is equal to that of the separate springs S-]_ or Sg. 

(5) Problems where component cases are coupled: The statements 

tinder item (4) are true only if the deformations for cases 1 and 2 can 
occur independently without restraining each other, such as the bending 
and shear deformation of a sandwich column with membrane faces. If this 
is not so, for example, if the faces have their own proper flexural 
rigidity, the internal action is increased above that from case 1 or 
case 2 alone. This can be taken into account in several ways. If the 
coupling occurs externally by a third agency, such as in the above example, 
this coupling case (case 0) can be first ignored; that is, the proper 
flexural rigidity of the faces is assumed to be aero. Finally the axial 
restraint offered by this coupling rigidity is added to the buckling load 
of the remain i ng system. This was done, for example, in reference 11. 

If cases 1 and 2 themselves are coupled, that is, if they are coupled 
internally (so that a deformation from case 2 adds to the internal 
reaction from case 1 and vice versa) as occurs, for example, in refer- 
ence 13, the extra restraint offered by case 2 is simply added to the 
internal action from case 1 and vice versa. Also a combination of these 
methods can be used. 

(6) Expression of internal actions and hence of buckling load in 

terms of buck.li.Tig loads for component cases: The internal action for 

case 1 due to a deflection w^ is equal to the external action for that 

case, if it occurs separately, which is proportional to the bu ckling 
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load P^ for that case and can he expressed In terms of P^ and w-j_. 

For example, for case 1 of the sandwich column considered in the section 
"Introduction," the internal moment from case 1 is equal to the 

external moment P^w-^. According to item (4) this is also the internal 
moment Ik for the combined case. Hence from item ( 3 ) the equation of 

external and internal moment for the combined case leads to the equa- 
tion P cr + w 2 ) = P-jW^. Similarly, the internal moment for the com- 
bined case is equal to that for case 2, which is again equal to the 
external moment for that case, giving Mj_ = PgWg. Equating this internal 

moment to the external moment gives P cr (w-j_ + Wg^} = P2 W 2* E limina tion 
of w^ and Wg gives P cr in terms of P^ and Pg. If the component 

cases are coupled and the latter method under item ( 5 ) is used, as will 
usually be done in problems where cases 1 and 2 have different boundary 
conditions, the internal action will contain both partial deflections w^ 

and Wg. Further, in such problems the external action is usually 
expressed as P cr ( w q + <PWg) or P cr (wg + yw-jJ where cp and 7 differ 
from 1 , so that after elimination of w^ and Wg a more complicated 
formula for P cr in terms of P^ and Pg results than for problems 
where cases 1 and 2 have the same boundary conditions. 

If splitting into three different cases is necessary the same method 
is used, equating the external actions to the internal reactions due to 
cases 1, 2, and 3 , respectively. This leads to three homogeneous linear 
equations in w 1 , Wg, and Vy Only when caseB 1, 2, and 3 have the 

same boundary conditions and are not coupled or are coupled by a fourth 
restraint that can be split off, as mentioned in the first case in 
item ( 5 ), does this lead to simple formulas for P cr as expressed in 

terms of P^, Pg, and P y such as those for composite columns coupled 

by elastic couplings with equal spacings (ref. 6). In other problems 
this method leads to a cubic equation for the buckling load P cr - How- 
ever, by combining the first two cases and then combining them with the 
third case (in ref. 13 such a case is presented, namely, the buckling 
of columns with T-sections), rather simple end formulas could be obtained 
containing square roots only. 

(7) Accuracy of method: The method leads to exact results when the 

deflections w n have the same shape, but sufficiently accurate results 

can be obtained if they differ in shape, as will be shown in this report. 

As shown in reference l4, if w-^j W£> . . . w q differ in shape, 
the method remains exact if P^, Pg, . . . P n are considered as the 
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buckling loads belonging to the shapes of w^ and Wg in the component 

cases, as they occur in the combined case. Since, in general, these 
are not the opt imum shapes that lead to the minimum buckling loads P-j_ 

and Pg as they occur in separate cases 1 and 2, by considering P^ 

and Pg as the minimum buckling loads conservative results will be 

obtained. However, if the structure is externally redundant, such as 
in the case of sandwich plates, there is another influence that may tend 
to make the method unconservative . This is the fact that in the composite 
case the boundary and continuity conditions have to be satisfied by the 
total deflection w only and not by the component deflections 
w i, Wg, . . . w n separately. This causes a relaxation of restraints 

that tends to make P-^ and Pg smaller for the composite case than for 

the separate cases and thus to make the result unconservative. Hence in 
such a case the method may slightly overestimate the buckling load P cr 

for the combined case. 

(8) Problems where buckling deflection in one individual case is 

arbitrary: Let the case in which the buckling deflection is arbitrary 

be case 2. Then in the composite case the deflection Wg for that case 

will have the same shape as w-^ for case* 1. This may be understood as 
follows: If during the deformation for case 2 the shape of w-^ and 

hence the internal reaction for case 1 does not change, the internal 
reaction for the combined case will also be determined by the shape of 
w-^ (from items (4) and (6)). Hence, in order to make the external 
action equal the internal one at each point, the total deflec- 
tion w^ + Wg, and therefore Wg, must also have the same shape as w-^. 

In some cases, such as that of a sandwich column with asymmetric boundary 
conditions which will be considered in this report, the shape of w-^ for 

the combined case will differ from that for the individual case 1. How- 
ever, also then, in order to obtain equal external and internal moments, in 
the combined case Wg will have the same shape as w^, but the shapes of 

w-^ and Wg will differ from that of w-j_ in the separate case 1. Exact 
results can be obtained by calculating P^ and Pg for the buckling 
deflections as they occur in the composite case. 

(9) Columns or plates with initial deflections: The method may also 

be applied to columns or plates with initial deflections. Here one of 
the cases is the initial deflection with an individual buckling load 
equal to zero (refs. 1 and 15 ). 

(10) Reduction to case of columns or plates on elastic foundation: 

In several problems it is convenient to reduce the restraint offered in 
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the component cases to that given by an equivalent elastic foundation. 
The structure then can be considered to be without any flexural or 
shearing rigidity but to be laterally supported by springs arranged in 
series with as many component springs as the number of cases into which 
its behavior is split. For example, if there are two cases, the springs 
consist of individual springs Sp and Sp (fig. l) . 


PROBLEMS WHERE ALL COMPONENT CASES HAVE SAME BOUNDARY CONDITIONS 


General Formulas 

If all component cases have the same boundary conditions the 
result ing formula expressing the critical load in terms of the critical 
loads of the separate cases usually acquires the general form 

F cr - p O + p r - p 0 + ( p l _1 + + V* + • • •) _1 U> 

where P-^, Pg, . . . P , are the critical loads of the separate cases 
and Pq is the critical load that is due to the rigidity from the coupling 
case. For example, for a sandwich plate Pq is the critical load from 

the proper flexural rigidity of the faces alone. Of course, this does 
not apply if only one mode is considered, such aB in the case of homo- 
geneous plates under nonhomogeneous stress distribution and that of the 
ultimate load of compressed plates, which is dealt with subsequently. 

For plates where the rigidities in the X- and Y-directions differ 
a great deal, such as corrugated- core sandwich plates, it is possible 
and often necessary to allow for ratios Wg/Vp of the component deflec- 
tions which differ for the imaginary strips running in the X- and 
Y-directions of which the sandwich plate can be thought to -be composed. 

This leads to a formula of the form 


P 


cr 



+ P 

r 





4 - • • 



(2) 


Here the subscripts x and y refer to the buckling loads of the plate 
due to the rigidity of the X- and Y-strips, respectively, also taking 
account of the influence of the torsional moments acting on those strips. 
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This formula leads to accurate results for orthotropic sandwich plates 
and also to more accurate results than equation (l) for isotropic sand- 
wich plates with different houndary conditions in the X- and Y-directions , 
such as long plates that are clamped at the unloaded edges. 


Discussion and Derivation of Equation (l) 

Equation (l), where in some cases Pq is zero, was derived in refer- 
ences **-, 6, and 11 for columns as well as for plates by applying 
items (l) to (6) and (8). In the case of columns the actions and reac- 
tions considered were the external and internal moments, and in the case 
of plates the transverse deflecting and restraining forces acting upon 
an element were compared. As stated in reference 11, this formula leads 
to accurate results if w^, w 2 , . . . v and w have the same shape. 

However, even if w-^, W 2 , . . . w and w are of different shapes 

sufficiently accurate results are obtained. In reference l4 the same 
formula was derived for the case in which w^ and w 2 are different in 

shape by comparing the work done by deflecting and restraining forces. 

It was found that equation (l) is exact if in the combined case w^_ 
and w 2 have the same shape as they have in the component cases 1 and 2. 

But this is obviously so only if, in both the component and the combined 
cases, w^, w 2 , and w all have the same shape, since then, from 

items (5), (M, and (6), in the combined case the distribution of the 
deflecting and restraining forces will be similar. That is, if in the 
combined case the restraining forces from case 1 are given by ^ > p w '] jn ®(x,y), 

where $(x,y) is a definite function of x and y, for case 2 they will 
be PpWp m <&(x,y) while the total deflecting forces will be P qr w m $(x,y) . 

Here x and y are the coordinates in the plane of the plate, w 1tt1 

and w„ are the maximum deflections from cases 1 and 2, and v = w, + 

2m 7 m lm 2m 

is the total deflection. 

Perhaps a better insight into the physical meaning of the method of 
split rigidities may be obtained by deriving equation (l) by using 
item (10). Consider a long rectangular sandwich plate compressed in 
the X-direction with arbitrary boundary conditions at the loaded and 
unloaded edges (fig. 2). According to item (l) the buckling deflection w 
is split into case 1, from bending, with deflections w^ and case 2, 

from shear deformation, with deflections w 2 - However, as shown in the 
longitudinal section in figure 3, the shear deflection c ann ot occur 
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without "bending of the faces (the bending rigidity of the core is 
neglected) . This means that in the combined structure case 2 c ann ot 
occur without inducing a partial deformation according to case 1; in 
other words, cases 1 and 2 are coupled. According to item ( 5 ) one way 
to deal with this situation is to split off the coupling action. This 
means that in cases 1 and 2 the faces are assumed to have no proper 
flexural rigidity. Since, however, the buckling deflection of the two 
single faces requires a load equal to their, proper buckling load Pq, 

from item (6) their deflection in the combined case will generate an 
extra load Pq on the plate. Hence, this load Pq has to be added to 
that from the combined cases 1 and 2. The combination of cases 1 and 2 
is denoted as the reduced case. This procedure is exact if the buckling 
deflection of the single faces (case 0) is similar to the total buckling 
deflection w in the combined case. 

In case 1 (deformation by bending only) the deflecting force acting 
on a small element H dx dy, where H is the total plate thickness 
(fig. 3), is 


D 1 " " P 1 


8x^ 


dx dy 


(5) 


where P-j_ Is the buckling load per unit width for case 1, neglecting 

the proper flexural rigidity of the faces. Hence from item (6) the 
restraining force is 


R, = -Pn 


8 2 w 2 


dx dy 


(4) 


If the plate was supported by an elastic foundation with a foundation 
modulus C-j_ the transverse restraining force acting upon an ele- 
ment H dx dy would be C-jW-^ dx dy. Hence the restraint offered by the 

plate in case 1 is equivalent to that of an elastic foundation with a 
foundation modulus (spring constant) 



( 5 ) 


In general, this spring constant may be a function of x and y. 
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Similarly, the equivalent spring constant for the restraint offered 
by the plate in case 2 (shear deformation of the core only) is 

C 2 ‘ - P 2 (6) 

where P 2 is the buckling load per unit width for case 2. Hence, if 
with the actual buckling deflection w the partial deflections w^ 

and W 2 have the same shape as in the individual cases 1 and 2, the 

reduced plate (without proper flexural rigidity of the faces) offers a 
restraint against buckling that is equivalent to the lateral restraint 
from two spring systems built in series with spring constants C-j_ 
and C 2 (fig* 4( a ) ) . Therefore, a unit lateral load per unit surface 
causes a deflection 



so that the equivalent spring constant for the reduced plate is 


C 


r 


C 1 C 2 _ 


Co + Cc 



( 8 ) 


The equivalent spring constant for case 0 (the single faces), with total 
buckling load Pq, is in the same way 


C 


0 



(9) 


Hence, if the actual buckling deflection w has the same shape as Wq 
( for case 0 alone), the equivalent spring constant for the actual plate 
is 


C = C, 


C r " 


Co + 


c C 2 

cTTc: = c o + 



(10) 


that is, the spring in figure 4(b), with spring constant C, is equivalent 
to the spring system in figure 4(a). 
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A plate with a spring constant C will buckle if for any element 

the deflecting force -P ^ is equal to the restraining force Cw, so 

3 xr 

that the buckling load is 




( 11 ) 


or, from equations (5), (6), (9), (10), and (ll), 

r 

'0/ 


A 


cr o 

8 2 


w J 


a 2 w 0 /ox £ 


0 


w. 


0 


P 1 P 2 (^l/^) { h \/ bx2 )/{ 

[- 1-2 

) 

p i(8 2 wi/8x 2 )/wi 

+ 

P 2^ 2w 2 j 

!^)i 

1-2 




0 aVa* 2 - w o 


P P 
12 


S 2 w/Sx 2 w 2 , r, d 2 w/dx 


( 12 ) 


3X ^ , p 


W 1 


With the aforementioned assumption that w = Wq, the first term to the 
right becomes Pq. The spring forces in the springs with constants 
and C 2 built in series are 

C^wq = CgWg (l?) 


so that from equations (5) and (6) 


Further, 


= P 2 

8 2 w 2 /8x 2 ^l 


W = W^ + Wg 




(15) 
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Using equations ( lA) and ( 15 ) , in equation (12) the denominator of the 
second term to the right becomes 


P 


1 


^2 

W 


P 2 


b%kx 2 V 1 

s\/a x 2 v 



^2 

v 


+ P, 



w 


- ( P 2 + P l) 


^S- 1 - V 1 

W 


= Pi + P 2 


(16) 


Hence equation (12) transforms to 


P_„ = P n + — 
cr Op 


P P 
12 

1 + P 2 


( 17 ) 


which is identical to equation (l) with P^ and so forth equal to zero. 

The more general equation (l) can be derived in the same way. Hence 
equation (l) is exact if w-j_ and w 2 differ in shape, but if in the 

composite case Wq = w, w^ and Wg have the same shape as in the 

separate cases 0, 1, and 2, respectively. This confirms the derivation 
in reference l4. However, as stated previously, this condition can be 
met only if Wq, w^_, and w 2 have the same shape. In the latter case 

equation (l) is exact within the limitations of the usual assumptions, 
such an neglect of the deformations from transverse shear in the faces 
and from normal stresses a z in planes parallel to the middle plane of 

the sandwich plate. 


Sandwich Columns With Various Boundary Conditions 

Calculations of buckling loads of sandwich plates with several 
boundary conditions were given in references 11, 12, and lk. In sand- 
wich columns or wide plates compressed in the short direction, the author 
considered until now only the simply supported case (fig. 5(a)). Here 
the deflection w^ for case 1, deformation from bending only, is a half 
sine wave. The buckling deflection W 2 from shear deformation alone 
(case 2) is arbitrary. Hence, from item (8), in the combined case w^ 
and W 2 will both have the shape of a half sine wave. The single faces 
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(case 0) will also buckle in a half sine wave, so that, from item (6), 
in the combined case the internal moment is increased by a sinusoidal 
moment PqW, where P Q is the buckling load of the faces. By adding a 

load P q to the buckling load P r of the reduced case an external 
moment PqW of the same amount is added, so that with a sinusoidal 
deflection w equilibrium is maintained in all cross sections. Hence 
equation (l) is exact for this case. 

In reference 23, the author's method was extrapolated for calcu- 
lating the buckling load of sandwich columns that are clamped at both 
ends or at one end (figs. 5(b) and 5(c))* It was found that for a strut 
that is clamped at both ends equation (l) yields exact results, but that 
for the simply supported clamped case (fig. 5(c)) equation (l) gives 

results that for two examples were unconservative by 5 and 6- percent. 

This is understandable from the foregoing discussion. In the 
clamped-clamped case (fig. 5(b)) the buckling deflections Wq and w-j_ 

for the individual cases 0 and 1 are both full sine waves and the shape 
of W£ is arbitrary. Hence, from item (8) and the above discussion, 
in the composite case w, w-j_ and W 2 will also forma full sine waves 

so that from item (l) the present method leads to exact results. Another 
way to see that equation (l) is exact for a clamped-clamped strut is to 
observe that the inflection points B and C are (l/4)Z and (3/4 ) l from 
the left end. The line of action of the compressive force P cr passes 

through these inflection points. The center part BC is in exactly the 
same condition as a simply supported strut with length L = ( l/2) Z and 
the left and right parts AB and CD can be considered as the right and 
left parts of such a column, respectively (fig. 5(b)). Since P cr , Pq, 

P-^, and P£ in equation (l) also apply to this simply supported column 

with length L, this equation is also exact for the clamped-clamped 
column. 

Consider now the clamped s imply supported column of figure 5(c). 

For case 1 the buckling deflection is part of a sine wave, with the line 
of action of the conpressive force as the axis. If the external and 
Internal moments Mg and are expressed in terms of the distance of 

the column to that axis, the equation Mg = Mj_ becomes (fig. 6(a)) 

P lV l = -Eglsd^/dX 2 (18) 


This equation and the boundary conditions are satisfied by 



NACA TN ^K)85 


19 


sin J x ( 19 ) 

where L is the effective length AB^ of the column (fig. 6(a)). The 
same applies to Vq for case 0. Since for case 2 the deflection w 2 

is arbitrary, from item (8) in the composite case it would have the same 
shape as w^_ and w. However, it may be easily seen that in the com- 
posite case w-^ has not the same shape as in the separate case 1. It 
is well known that in asymmetrically restrained, transversely loaded 
spans the distribution of the bending moments is changed by shear deflec- 
tions as, for example, occurs in thin tubes or shells supported at three 
or more points. The effect of clamping a support diminishes . From fig- 
ure 6(b), presenting the composite case, it is seen that the line of 
action AC 2 of the compressive force embraces an angle a with the original 
axis AC of the sandwich column. Hence at the clamped end C the trans- 
verse shear and the shear deformation is not zero, so that, although in 
case 1 (fig. 6(a)) the slope at C is zero, in the composite case the 
slope angle 0 at C is not zero. Now, considering first the reduced 
case (P 0 = Oj, from the same reasoning as given in item (8), w-j_, W 2 , 

and w will have similar shapes with respect to the line of action of 
the compressive force but different shapes with respect to w^_ in case 1. 
Hence, from items (l) and ( 7 ) and the foregoing discussion, equation (l) 
is not exact if ^ and P 2 are considered as the buckling loads for 
the separate cases 1 and 2. On the other hand it will be exact if P-^ 

and P 2 are the buckling loads for cases 1 and 2 as they occur in the 
reduced case of figure 6(b ) . 

In the same way as shown for case 1 in figure 6(a), in figure 6(b) 
the deflections w^ for case 1 are given by equation ( 19 ), so that w 2 

and w will also vary sinusoidally with respect to the lin e of action 
ABC 2 of the compressive force P cr . Since at C the elastic line of the 
strut as shown by the solid curve embraces an angle with its original 
axis AC, the effective length L for the reduced case (fig. 6(b)) is 
greater, with respect to l , than in the separate case 1 (fig. 6(a)), 
where, as is well known, L = 0.71* Therefore, In the actual case P-j_ 

is smaller than it is for case 1 in figure 6(a). Since P 2 Is inde- 
pendent of shape or length, it is the same as in the individual case 2. 

In reference 25, P]_ from the separate case 1 (fig. 6(a)) was used in 

equation (l), which explains why P cr was found to be higher t han the 
exact value (Pq was assumed to be zero^ . 

The effective length in figure 6(b) may be found in a similar way, 
as was shown on pages ’JO to J8 of reference 6 for various problems, 
using the Haarman method. The deflections w-^, w 2 , and w are measured 
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from the lin e of compression ABCg- With the deflection at the 

clamped edge no rotation can occur with respect to the original axis AC, 
so that the tangent C^'C^ is parallel to AC. Hence from figure 6(h) 


a = 




( 20 ) 


where w^ is given by equation (19) and, since v has the same shape 
as was measured from ABC2> 


w = v m sin | x 

Insertion of equations (19) and (2l) into equation (20) gives 


( 21 ) 


( 22 ) 

L Win L 

From items ( 3 ), (4), and (6) the equality of external and Internal 
moments requires that 

P r w = P^ = P 2 w 2 (23) 


so that 


= Zi = y l = P 2 

W P-j^ + Pg 


(24) 


Since P-^ is equal to the P^ for a simply supported strut of length L, 
from reference 11 



(25) 


where (fig. 3) 

I s = ( 1/2 )h( t + h) 2 


(26) 
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and 


p 2 


(t + h) 2 
t G -c 


(27) 


where G c is the modulus of rigidity of the core. Hence, from equa- 
tions (22), (2k), and ( 25 ), L has to he calculated from the following 
equation: 


* l 2 

tan jt T = -= 

L it?ELI 


rtP^Z/L 


V 



1 1 + P. 


( 28 ) 


Then P-^ and P 2 follow from equations ( 25 ) and ( 27 ) and P cr = P r from 
equation (l), with Pq = 0 and P^, . . . P Q = 0. With the dimensions 
used in an example in reference 23, E S I S = ^3*200 kg-cm 2 , P^ = 360 kilo- 
grams, and l 2 = 800 cm 2 , equation ( 28 ) yields l/L = 1 . 276 , so that 
L 2 = i^/l .628 = 491 cm 2 . Hence equation (25) gives Pj^ = 870 kilograms, 
so that from equation (l), with Pq = 0, P cr = 25^ kilograms. This is 

in accordance with the result obtained in reference 23 from an exact 
calculation, which shows that also for columns with asymmetric edge condi- 
tions equation (l) is exact if applied to the correct effective length of 
the column. However, for such cases, it loses the advantage of its 
simplicity. 

Figure 6(h) illustrates what was said in reference 14- and under 
item ( 7 ) • The deflections w^_ and Wg satisfy the boundary conditions 
together and, as seen clearly from this figure, w^ does not satisfy 
the boundary conditions . This relaxation of restraints causes P-^ to 
be smaller than P^ from the separate case 1 alone (fig. 6(a)). 

Usually Pq is negligible, but if the bending rigidity of the faces 
is taken into account, one simply can calculate Pq in equation (l) for 
the effective length L from equation ( 28 ), so that from reference 11 


P 


0 


rtSlgh 5 

6 (l - v 2 )l 2 


(29) 
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Actually this makes equation (l) slightly conservative, since from fig- 
ure 6(b) the deflection w has a sharp break at the clamped edge, so 
that Fq actually will be slightly higher than would follow from, 
equation ( 29 ) . 


Columns With Batten Plates and Latticed Columns 

An application of equation (l).that shows some other features is 
the calciilation of the critical stress of metal columns connected by 
batten plates with equal spacings (fig. 7(a)) or of timber columns that 
are coupled at equal distances (refs. 3 and 6). In many applications 
the deformation of the batten plates can be neglected. Then the buckling 
deflection can be split into two cases. Case 1 is caused by bending with 
respect to the common axis (fig. 7(^})- Case 2 is the deformation by 
shear, which bends the single columns between the batten plates in S-curves 
(fig. 7(c)). These deformations can occur independently of each other, so 
that 


P 


cr 




-1 


(30) 


fcy 


Including the plastic range, 


P-j_ is sufficiently accurately given 


Pi = 



(31) 


where 


I = 2I c +(Ah 2 A) ' (32) 

Here is the tangent modulus, h is the spacing of the axes of 

inertia of the single columns, A is the total cross section of the 
composite column, and I c is the moment of inertia of the single col umns ■ 
In case 2 the batten plates translate with respect to each other without 
rotation (fig. 7(c)), so that 


o t 



(33) 


where c Q is the effective free length of the single columns between 
two batten plates. 
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If also the deformation by shear and bending of the batten plates 
or wooden connections has to be taken into account* case 3 is added 
(fig. 7(<i))- However, now cases 1 and 3 cannot occur independently, 
since the deflection w^ for case 3 bends the single columns in a single 

half wave with respect to their original axes, which is part of the defor- 
mation in case 1. Therefore, according to item (5), the proper rigidity 
of the single columns is first assumed to be zero and subsequently its 
influence is added in the form of Pq. Hence, in equation (3l) for Pq 

the moment of inertia I is that of the reduced col umn . 


Ir = Ah 2 / 4 (34) 

since I c is assumed to be zero) P 2 remains as given by equation (33) . 
From reference 6, assuming only the batten plates (or wooden connections) 
to deform, the individual buckling load for case 3 is 

P 3 = 0h 2 /cv ( 35 ) 

where c is the center-to-center spacing of the batten plates and v is 
the translation of the single struts with respect to each other per -unit 
shear force Q (fig. 7(d)). In reference 6 the factor 0 was calcu- 
lated for several cases. If the rigidity of the end batten plates is 
half that of the intermediate ones, 0=1. 

In the present problem Pq is not equal to the buckling load of 

the single struts (in contrast to the case of a sandwich column, where 
it was equal to the buckling load of the faces), because during the 
deflection W 2 from case 2 (fig. 7(c)) no deflection in a single half 
wave occurs. The bending in S-curves between batten plates causes many 
inflection points to occur in the single struts, so that the internal 
moment in them is not increased. If the half-wave bending of the single 
struts occurred simultaneously with its total deflection w, from item ( 5 ) 
it would add a resistance in the longitudinal direction of the strut of 

2jt%1(-I g /z 2 . However, since this half-wave bending occurs during the 
deflection w - Wg only, its curvature reduces in the ratio - wjj w, 

so that the added axial resistance is 


w 


P 0 = 


- w 2 2x%.l c 


w 


(56) 
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Similarly to equation (23) 


P 2 W 2 = p r w 


(37) 


so that 



(38) 


and 



(.39) 


Since with a sufficiently large number of "batten plates w^ is distrib- 
uted sinusoidally and since for w 2 and w^ the buckling deflection 

at the batten plates is arbitrary, by the same reasoning as given in 
item (8) for one arbitrary component deflection in the composite case w, 
w 2 , and Wj will also have a generally sinusoidal half-wave distribution 

(actually, from ref. 24, if c Q = c, the centers of the batten plates are 
situated on a half sine wave). 


From item (10) the column can be imagined as being without flexural 
rigidity and being laterally supported by a spring system as shown in 
figure 8. The equivalent spring constants are given by equations (.5), 
(6), and (9) and by a similar equation for Cj. In the same way as shown 

for a spring system with spring constants Cq, C^, and C 2 , this leads 
to an equation like equation (l) which now contains a term with P^. The 

requirements of items- (l) and (7) for equation (l) to be exact are not 
entirely fulfilled here, since between the batten plates the elastic 
lines in the various cases are not exactly similar, although their over- 
all shapes are those of half sine waves. They will be closer to that 
shape, the larger the number of batten plates. Equation (l) for this 
problem is 


r- 


P cr P 0 + P cr P 0 + 



(4o) 
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where 


Pq, T-i, P 5 , and P^ are given by equations (39)> (3l) with. I 


■ 2 > 


replaced hy I r from equation (34), (33)> and (35)* Dividing by the 


cross section A gives 


<?cr _ °b + a r - + 


A 


= *0 


(° 


_l\ -1 

+ [a-, - 1 - + a 2 + P3 J 


-1 ^ „ -1 


° ilr/rf 

where l r j r is the effective slenderness of the reduced case. 


jt 2 E, 


ff-, = 


1 (A 


A 


ff 2 \2 


( C o/ r e)' 


a - pfa2 
3 cvA 


(4l) 


(42) 


and 


r 2 = I r /A = h 2 /4 

From equations (39 )j (4l), and (42), 


„ _ A M * E t 
0 \ ^ n2 


(W 


1 - 


/ c o/ r c 

,V r 


A 




(43) 


(44) 


As shown more extensively in reference 6, equation (l(-l) may be written 
as 



(45) 
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where the effective slenderness Z e q_/ r given by 


M 2 . 

( l r/> 

^) 2 

r r 

N ' 1 + 

( 2 r/ r ) 2 


'o/*c) 2 ( T c/ 1 ) 2 


(46) 


From equations (4l) and (42) 


( 2 r/ r ) 2 = (Vr) 2 + 



it 2 E.j.Acv 


/K) 


(W 


With connections of vanishing rigidity, the factor v in equation (47) 
becomes infinite, so that equation (46) yields the correct result, 
namely, that the effective slenderness is equal to that of the single 
struts l/r c . In contrast, for example, equations (348) to (350) of 
reference 25 and equations (96) and (97) of reference 26 for infinitely 
weak connections yield zero buckling loads. From figures 7(c) and 7(d) 
a batten plate at the center of a col umn has no effect, so that for a 
column with end connections only, one has to assume that the spacing of 
the batten plates is c = Z/2. 

An exact formula for the buckling stress of these columns has been 
derived in reference 24 for the case in which c 0 = c and the rigidity 
of the end plates is half that of the intermediate ones. Changing the 
reduced modulus to the tangent modulus Hj. gives the buckling condition 



where n = Z/c and m = (^eq/ r ) /( c / r c)* 4his same case, in equa- 

tions (46) and (47) c c = c and 3=1. Table V of reference 6, given 
as table 1 herein, contains some values of m calculated from equa- 
tions (46) and (47) as well as from equation (48) by Mr. Lie Han Yang. 

It refers to timber columns built up from single struts with dimensions 
8 by 20 centimeters, the smaller dimension being that parallel to the 

plane of bending; E t wan assumed to be 100,000 kg/cm 2 . Table 1 shows 
that equation (46) is conservative and indeed is more accurate for large 
numbers n than for small ones. The general accuracy is very satisfactory. 
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Similar formulas were derived in reference 6 for latticed columns, 
coupled by diagonals only (fig. 9(a)) or by diagonals and verticals 
(fig. 9(4*)), using equation (l) with P^, . . . P n equal to zero. If 

I is replaced by I r from equation (34), Pq is given by equation (3l)i 
Pq represents the axial resistance of the single col umns , so that 


P 0 - 


2jt' 




(49) 


In order to find Pp it is observed that for case 2, where the single 

columns are assumed to be infinitely rigid against axial strain, a 

slope 0 of the column with respect to its original axis causes a tr ails - 

verse shear force Q = P 2 P that has to be resisted by the lacing . Denoting 

as Q' the fictitious transverse shear force that would cause a unit 
angular distortion, the equation Q = Q'0 is obtained, so that 


P 2 = Q* 

Inserting P Q , Pq, and Pg into equation (l) gives 


(50) 




p cr ~ o 


2I P + 


+ 


(m2E t I r /Q') 


(51) 


As may be easily checked, for a column with diagonals only (fig. 9(a)) 


Q’ = EAq sin^a cos a (52) 

where A<q is the cross section of the diagonals. For the arrangement 
shown in figure 9(t) 


1/Q’ = 


-/( E Ad 


. p 

sxn^a 


cos a 




^tan 


a/EA v j 


(53) 


where Ay. is the cross section of the verticals. For both these cases 
an upper limi t for P cr is that of the single struts with a free 
length c and slenderness c/r s . With diagonal cross sections Aq = 0, 
equations (52) and (53) yield Q’ = 0 and equation (51) reduces correctly 
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to the critical load of the single columns, in contrast with equa- 
tions (335) to (339) of reference 25 and equations (93) and (9^) of 
reference 26 which yield P cr = 0 . 


Sandwich Plates With Orthotropic Core 

A comprehensive discussion on the application of equation (l) to 
the elastic and plastic buckling of sandwich plates was given in refer- 
ences 11 , 12 , and 1 ^. It was shown that equation (l) is exact if Wp 

w 2 , and w have the same shape, such as occurs for a compressed long 

plate with simply supported unloaded edges. In that case the require- 
ments for exactness, stated in item (l), are satisfied. For a long plate 
with simply supported long edges that is subjected to shear, the deflec- 
tions w-^ and w 2 differ somewhat, but nevertheless results obtained 

in references 22, 23, and 27 differ less than 1 percent from those from 
equation (l). This is due to the fact that, although w-j_ and w 2 

differ in shape, the rigidities of imaginary X- and Y-strips (parallel 
and perpendicular to the long edges) are still of the same order of 
magnitude "being simply supported at the nodal lines and at the long edges, 
respectively. 

A more i nt ricate case arises if for one of the component cases the 
strips manning in the X-direction (the X-strips) have a rigidity that 
differs from that of the Y-strips (fig. 10 (a;)). For isotropic sandwich 
plates this occurs for a long compressed plate that has clamped unloaded 
edges. For case 1 (deformation from bending) the clamped Y-strips are 
much more rigid than the X-strips, which can be assumed to be simply sup- 
ported at the transverse nodal lines. On the other hand, in case 2 
(deformation by transverse shear forces) in any buckle both X- and Y-strips 
deflect in half sine waves so that their rigidities are of the same order 
of magnitude. As explained in reference l 4 and item ( 7 )> for cases where 
w-]_ and w 2 differ in shape two influences occur, one that tends to make 

equation (l) conservative and another that tends to make it unconserva- 
tive. In the present case the latter influence, due to the fact that 
the boundary and continuity conditions have to be satisfied by the total 
deflection w = w-^ + w 2 and not by w-j_ and w 2 separately, is pre- 
dominant, making equation (l) unconservative by 7 percent (ref. 28). 

Apparently the fact that in the combined case w^_ and W2 do not 
satisfy the boundary conditions separately is mainly due to the differ- 
ence in rigidity of X- and Y-strips and not to the difference in shape 
of the deflections w^ and w 2 . By prescribing a constant ratio ^m/^lm 

for the amplitudes of the deflections in cases 2 and 1, one prescribes 
that the ratio w 2 xm/ w lxm am Pli'k U( 3 - eB of the center X-strip'is 
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equal to ttie ratio ^2ymj w lxia. o£ " fciie Moplitudes of the center Y- strip. 
Actually in the composite case this cannot he true. If w 2xm/ w lxffi. ds 
equal to cl, w 2ym/ w lym ^^1 be much more than a because of the 
clamping of the Y-strips since w-j^ is and Wg^ is not affected by 
the clamping, as explained’ previously. Hence, with a prescribed ratio 
W 2m/ W M = w 2xm/ w lxm = w 2ym/ w :Lym> one could iJ^ag^e fo ? tlie center 

Y-strip Wg-y. has negative end deflections ~ w 2p v ly das positive 

end deflections (fig. 10(b)). Then the boundary conditions 

at y = 0 and y = b are satisfied by w = w^_ + = w-^ - w^ = 0, 

while the ratio w 2ya/ v lya amplitudes of the actual deflections 

is more than the ratio v 2ya/ w lym f as ii should be. This shows the 

relaxation of restraints for the component cases 1 and 2, which makes 
the actual buckling load smaller than that obtained from equation (l). 

Although for isotropic plates with different rigidities of X- and 
Y-strips equation (l) is still sufficiently accurate, for sandwich plates 
with anisotropic core, where the ratio Gyj G x differs too much from 1, 

as stated in reference l4, equation (l) becomes too inaccurate. From 
the preceding discussion, the obvious way to improve this situation is 
to admit a ratio w 2ym/ w lym i° r ’ fcde Y-strips that differs from the ratio 

w 2xm/ w lxm ^ or X- strips . This leads to equation (2) that will be 

shown to give accurate results for anisotropic sandwich plates and also 
more accurate results than equation (l) for isotropic sandwich plates. 
Equation (2) will now be derived. 


Derivation of Equation ( 2 ) 

Consider a long sandwich plate with isotropic faces, an orthotropic 
core compressed in the long (X) direction, and arbitrary boundary condi- 
tions at the unloaded edges (fig. ll). After splitting off the proper 
rigidity of the faces (item ( 5 )) the plate is considered as a grid 
consisting of X- and Y-strips, but due account is taken of the restraints 
exerted upon these strips by the twisting moments along their long edges. 
This same method was applied by the author in reference 29 and in earlier 
papers for deriving simple formulas for the bending moments in rectangular 
plates under several loading and boundary conditions. From item (l) the 
deflections of these strips are split into those from bending about the 
common middle plane (case l) and from transverse shear (case 2). 

Let the lateral restraint offered by the X- strips in cases 1 and 2 
(cases lx and 2x) and that of the Y-strips in cases 1 and 2 (cases ly 
and 2y) be equivalent to that of elastic foundations with foundation 
moduli of C^, Cgx, C^, and. Cgy, respectively. Hence, if in the 



30 


KACA TN 4085 


composite case the deflections from cases 1 and 2 have the same shape as 
in the individual cases, from equation (8) the equivalent spring constant 
of the reduced plate is 


c r = o x + c y = ( 0lx _1 + o^- 1 )' 1 + (Oj y 1 + V 1 )' 1 (54) 

because the resistances of X- and Y-strips just act additively. The 
equivalent spring constant of the actual plate which is similar to 
equation ( 10 } is 

C = C Q + C r = C Q + ^C^" 1 + Cg/ 1 ) 1 + (Ciy" 1 + Cgy" 1 )' 1 (55) 

where C Q is the equivalent spring constant, for the proper rigidity of 

the faces (case 0). The equivalent spring system is shown schematically 
in figure 12. Using equations similar to equations (5), (6), (9), (13), 
and ( l4) , where the subscripts 1 or 2 are replaced by lx and ly or 2x 
and 2y, in a way similar to that in which equation (10) reduced to equa- 
tion (1), equation (55) reduces to equation (2). In equation (2) P Q , 

? 2 x’ P ly' 841(1 P 2y are the "buckling loads for the individual 

cases 0, lx, 2x, ly, and 2y, respectively. Similarly to equation (l), 
equation (2) is exact if the component deflections w 0 , w 1y , w^, w-^, 

and W2y for these cases have the same shape and if these shapes are 

the same in the composite case as in the separate cases. Since for 
case 2, for simply supported as well as for clamped edges, the X- and 
Y-strips deflect in half sine waves (refs. 5 and 12), equation (2) will 
be exact for simply supported isotropic or orthotropic sandwich plates, 
where for case 1 the deflections have this same shape. However, as will 
be shown, for other boundary conditions also equation (2) is very accurate 
(maximum discrepancy about 3 percent). 


Application of Equation (2) to Orthotropic Sandwich Plates 
With Simply Supported Edges 

Equation (2) can be applied to orthotropic sandwich plates with 
simply supported edges as follows. The restraining force exerted by an 
X-strip per unit plate surface can be written as (fig. 13 ) 
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where Q x is the transverse shear per unit width of the strip. From 
equation (102) of reference 30 , for case 1 ("bending only ) , by adding the 
subscript 1 , 




^ ^ 
s "r 


lx Sx Sy 

and in the same way from equation (37) of reference 30 


(57) 


(b%. 


“lx - 


-N c 


lx 


a 2 w. 


ax 2 


+ V 




ay 2 


(58) 


where N s is the flexural rigidity of the reduced sandwich plate, and v 

is Poisson’s ratio. As explained previously for the reduced case the 
proper flexural rigidity of the faces is neglected. The bending rigidity 
of the core is always negligible, so that (fig. 3) 


N s = 


E s h(t + h) 2 
2(l - V 2 ) 


(59) 


where E s is the modulus of elasticity of the faces ( skin) . From the 
equations on page 11 of reference 31 it follows that, using again the 
subscript 1 for case 1, 


M. = _m, « -35, ( ■&_ _ 1 8zx _ 1 azy) 

n lyx w Ixy ^Ixy ^ x £y 2 ay 2 hx.J 

The shear angles yx and yy can be expressed in terms of the deflec- 
tions Wg^. and Wgy for case 2 occurring from shear alone, since 

obviously 
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Further, since the faces are isotropic, the ratio of the twisting 

moment to twist is N s (l - v) as it is for an isotropic plate, so that 


M- 


lyx 


= -NJ 1 - v) 


8 2 w 


1 3 w 2x 


d 2 w, 


2y 


3x by 2 bx by 2 bx by/ 


(61) 


Since 


v 2x = 


= w - w- 


lx 


v 2y ~ w ‘ w 2y 


(62) 


M. 


lyx 


fb% b 2 w 


(63) 


Insertion of equations (58) and (63) into equation (57) and- of the latter 
into equation (56) yields for case 1 


R 


lx 


a \vx 

8x 2 ^T8y 


N c 


>4 
0 w- 


n4 

O W- 


hx 


+ i (1 _ V ) 

bx^ ^ bx^by^ 


1 b \r 

+ I (1 + v) — 

^ 8x^dy^ 


(64) 

If the lateral restraint given by the X-strips in case 1 was given by an 
imaginary elastic foundation with a foundation modulus C^, R^x would 

be equal to Cix w lx so that the equivalent spring constant is 


'lx 


'iy 


R la 

W s 

B Vx 

w lbi 

w ii 

8x^ 

the equivalent 1 

w 

£ 


.4 

a w, 
iy 

w iy 

W 0y 



4 4 

, 8 w 1v b ¥, 

+ I (l _ v) — + - (l + v) ^L. 

2 ^. 2^2 2 ^. 2>.2 


8x £ 8y" 


S^w, 

iy 

bx 2 b/ 


bxay 11 


5^w, 

lx 

8x 2 8y c 


(65) 


( 66 ) 
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The moduli of rigidity of the core against transverse shear forces 
Q x and Qy are denoted as G x and Gy, respectively. From equa- 
tions (6) and (7) of part I of reference Qx and Qy for case 2 
can he expressed as 


^2x 


*2y 


ft + h ) 2 p ^ 2 x 

” t x 8 x 

= ^ + h ) 2 Gy ^ 

t y 8 y 


■\ 


( 67 ) 


J 


so that, using equation (5 6) and analogous to equations (65) and ( 66 ), 

■\ 


"2x 


_ R 2x _ 1 5( ^2 x _ (t + h) 2 ^ V ; 

v 2x ~ v 2x ^x ‘ t V2x Sx 2 


2 x 


c. = ^ 1. = - ( * + *l f i 

w 2y v 2y * w 2y 


( 68 ) 


J 


For a long plate with simply supported unloaded edges the boundary 
conditions are satisfied if 


v w 0 w lx w 2x v ly ^y , it , it 
— = = — = — -i- = gin _ x sin t- y 

w Gm w lxm w 2xm v Iym w 2ym a ® 

Inserting this into equations (65), (66), and ( 68 ) gives 




■ 5 t K * 


°2y " T B = 
b 


-x + i (1 - v) -i- + J (1 + v) 


,4 2 


2 2 


9 xP £ 


1 + ^( 1 - v) 4 + 1(1+ v) 


2 2 


9yP 




(69) 


( 70 ) 
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■\ 


02x . £ L±±*r. Sx a 

^ t 2 t X p 2 


Co,, = 4 ■ (t - ?- ) .- G, 


' 2y ” t 5 t 


( 71 ) 


where a is the half wave length in the X- direction, 

P = a/b 

and 


■\ 


w- 


e 


* a T! 


Ixm 


m 


w- 


lyrn 


w. 


m 




(72) 


(73) 


From figure 12, since the springs with constants and Cg^ have 

to carry the same load, 


C lx w lxm “ C 2x w 2xm 


(74) 


so that 


w 


e v = 


lxm 


'2x 


X w lxm + v 2xm C lx + C 2x 


(73) 


Similarly 


w. 


0y - 

« v, + v, 




lym " r "2ym C ly + C 2y 


(76) 


Inserting and Cg x from equations (70) and (71) into equation (75) 

gives 


e„ = 


1 + 



1 + v 


(77) 
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where 

jr^ot 

r £ 

X + h ) 2 

Equation (77) can be written as 



From equations (7O), (71), and (76) 



where 

rt^TSLt 

p = - 

^ b 2 Gy(t + h ) 2 

Equation ( 80) is written as 


(78) 


(79) 


(80) 


(81) 
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In equation (2) P Q is the huckling load of the faces if buckling alone 
occurs, so that, from equation (g) on page 329 of reference 26 


. _ 2it^N 

0 = b 2 



m 


where N is the flexural rigidity of each of the faces. Further, in 
equation (2), in connection with equations (5), (6), (69), and (75), 



P P 
r lx 2x 

p lx + p 2x 


'2x 


c lx + C 2x 


P lx 9 X P lx 


(85) 


Similarly 

(*: ly ' 1 + V 1 )' 1 = W 

From equations (5) and (69) 




P lx = G lx 


p ly - ^ c ly 


Hence equation (2) can be written as 


p - 2m N 
P cr b 2 


p + P ) 2 + ^ ( 6xClx + 6yCly ) 


( 86 ) 


(87) 


( 88 ) 


which result could also have been obtained from equation (35) in combination 
with equations (9), (ll), (69), (75), (76), and (84). Using equation (70), 
equation (88) yields 


P 


cr 




(89) 
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where 0 X and 0y are given by equation ( 83 ). As stated previously, 
this result is exact. This will he shown in the next section by compari- 
son of a more general formula with the results obtained in reference 32 . 

For isotropic core G x = Gy = G and r x = ry = r, so that from equation ( 83 ) 


B 2 , 1 - v 
P + o 


©X = 9y = 9 = 


b 2 + 


(1 + P 2 )r 


(l + P 2 ^r + + p^ 2 


and equation ( 89 ) becomes 


(90) 


P cr - + A * eN s) (91) 

This result, can be shown to be identical to that given by equation (31) 
of reference 11 for t] = 1 (elastic range). 


Plastic Buckling of Corrugated Core Sandwich Plates, Hinged 
or Clamped at Unloaded Edges 

A cross section of the plate is given in figure l4. Considering 
first case 1 (bending only), from equation ( 3 °) of reference 33 or 
equation ( 22 ) of reference 34 the bending moment My from bending about 
the middle plane of the sandwich and carried by the faces, considered 
as membranes only, is 


^ly ”-®s^s 

where (fig. 3 ) 

Ib = § h (t + t ) 2 (93) 

For obtaining M^ the moment ) taken by the corruga- 

tion cross section, where E^ c is the tangent modulus of the corrugation 

material, must be added to the moment carried by the faces. Hence, using 
again the same equations from references 33 or 34, 


-s 2 

0 w- 


B 


lx 


S 2 w- 


dx 2 


+ D 


3yl 


ay 2 


( 92 ) 



38 


MCA TN 4085 


where 


S 2 w. 


M' 


•lx = 


" EsIs ' A Sx 2 


lx 


a 2 w. 




B 2 w. 


+ B - Efc c I c 


lx 


Bx 2 


- “Eglg 


/ E,\ B 2 w-, B 2 w, 

I A + p — 1 — + B 

l E c/ Bx 2 By 2 


P = 


El 
c c 

Eb^s 


(9^) 


(95) 


and E c is the elastic modulus of the corrugation material. From the 
same references the ratio N s (l - v) of the twisting moment to twist 
in the plastic range changes to 2E S I B F so that, instead of the present 
equation ( 63 ) , 


^lyx “ - E s I s F 


B 2 w n B 2 w.. 
lx . ly 


(96) 


Hie plastic parameters A, B, D, and F are given By equations (22) 
to (24) of reference 33 and equations (21) or (21a) of reference 34. 
Using equations (94) and (96) gives, instead of equation (64), 


P = . 

^ ^ x 2 Bx By 


= E s*s 


A + p 


Etc 


^4 

s w lx 


^4 

0 w. 


+ B 


..4 .4 


Ec J Bx^ Bx 2 By 2 \Bx 2 By 2 Bx 2 By 2 / 
from which, instead of equation (65), 


R 


■'lx 


w. 


Ox 

'lx 


Vs 


w- 


lx 


A + p 


Etc^ * w lx 


E, 


Bx* 1 


+ F 


B w- 


lx 


bV 


(97) 


Bx 2 By 2 


W-1, r 


> 


(98) 
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In the same way, using equations (92) and ( 96 ) gives 


f 

a 2 i O) 

/ El 

/w - 6 s 

" a 

O v ly 

T) -sy , -p 

sSf-, 

^ +(v, 4- y) 

^ w lx 

[ Sy 2 

Sx Sy J j 

r V " w^ 


Sx 2 ^ 

dx 2 Sy 2 


(99) 

Let the deflections now be, as was done for laterally loaded slabs in 
reference 35 and for a stability problem in reference 36 , 


where 


w^ = Y x sin | x 

. 

w- = Y v sin — x 
ly y a 

's 

*x=X w xmi u m 

► 

Y y = ^ ' w ym u m 

J 


( 100 ) 


( 101 ) 


Here and w^ are constants and the are the normal functions 

that automatically satisfy the boundary conditions at the unloaded edges, 
since they represent the modes of vibration of a Y-strip. These functions 
have the general form (ref. 37) 


Un = C! (cos ctjjjy + cosh o^y) + C 2 (cos - cosh cc^) + 


C 5 (sin cc^y + sinh + C^(sin ct^ - sinh a^y) (102) 
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so that 



>1 


- 


2 - 

— -— - ccjn. um 

dy 2 




(103) 


However, since the functions u m are orthogonal, ^ u^ can be 
developed in a series of Uj^: 


% * k lja u l + k 2m u 2 + k 3m u 3 + 


(104) 


It appears sufficiently accurate to use only the first terms of equa- 
tions (lOl) and (104). Hence 


^ U 1 2~ 2 
— 2~ = ^l^l = a l k H u l 
dy 


(105) 


Insertion of equations (100), (lOl), ( 103 ) , and (105) Into equation (98) 
gives 


E I c 

n _ S S 

c ^~~yT 


( A+ ^)s Y *- r 


2 d 2 ! 2 d 2 ! - 

— - (B + F) - 

g2 ,jy2 a 2 ^.2 


A 

- r E s l s 

b 4 6 6 


E. 


< A + d 


tc\ 1 

w 


b 2 a 2 k nn r 


"1 "11 


* 2 g 2 


9, 


F + (B + P) — 
s x 


J 


(106) 


Similarly one obtains from equation (99) 


C ly “ ”4 E b I s' 


foV - bV 2 k,„ r 


D 


1 11 
jt 2 B 2 




e 

F + (B + F) — 

e y 


■\ 


j 


(107) 
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Instead of calculating a,^ and. various bo undar y conditions , it 

may be observed that, if one or two of the unloaded edges are either 
simply supported or clamped, from equation ( ^fl) of reference 12, the 
buckling load of a homogeneous plate is 


p » • 7 EIi 




+ p(B + 2F) + qDfB* 1 


(108) 


where p and q are given in table 1 of reference 11. Hence, with 


w = Y sin — x 
a 


(109) 


where Y is a function of y alone, the equivalent spring constant is, 
analogous to equation (5)> 


n 2 / 2 k 

— — . O W" / jt -i— . Jt T-J -J- 


4- + p(B + 2 f)4 + ID 

P 


(no) 


From equations (106) and (107) for a homogeneous plate, where 0 y/ /0 x = 1 
and p=0, 


C - + C ly - ^ ®s^s 


JL 

P 11 - 


Hb^o^k. 


a. 


* 2 P 2 


(B + 2F) + 




(in) 


Comparing equations (no) and (in) gives 


*\ 


2b 




n 


= -P 








= q. 


(112) 
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so that equations (106) and (107) become 


"lx - 3 ■ c, s-‘-s 


C ly - ^ Vb 


r. \ 



■ L + n ! m ) 1 + JS- 

0 

F + (B + F) 


\ B c / 2p 2 

9 x 


V. 


y 


(113) 


r 


qD + 


F + (B + F) JE 


2£" 




For example., for clamped edges, from table 1 of reference 11, p = 2.5 
and q = 5. From references 35 and 3 6 for these boundary conditions, 
an = 2.365/c, where c = b/2 so that a^b = 4.73. Further, 
k n 1 = -0.54-984. With these values 




2b 2 o-, 2 k, , 

-L-L = _p c: 

*2 5 




= 5.15 


(n-4) 


which is in excellent agreement with the more accurate values -p = -2.5 
and q = 5 from equation (112). 

Defining, as in reference 31, the shear Btiffness Dq^ as the ratio 
of shear to shear angle or 


In case 2 


~ ^x/ 7 x 

D 0y = ^y/ 7 y 


8 w, 


■\ 


Q 2x “ D Q/x " D Q, x 

8w, 


2 x 


^ T> 


2y 


J 


(115) 


(116) 
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■where Dq^ and Dq^_ are given in reference 38 . This gives, analogous 
to equation (68), 


= 5&S = _1_ = ^ d^Wgx 

w 2x v 2x ~^ r w 2x 8x 2 

D 0y 8 2 wg^ 


Using equation (69) this hecomes 


C - ft 2 p 
c 2x a2 d Q x 


C 2y = 72 D 0y 


(117) 


(118) 


Insertion of the first equations in equations (113) and (ll8) into equa- 
tion (75) yields 


1 + r * A + 


-L. 

+ e[] 
c / p 2 2 L 


F + (B + F) tX 
°x 


(119) 


where 


*% I 

s s 




( 120 ) 


from which 


JC / 3 


1+ f A + \ + |F r x L0x + | (B + F)f x 9 v = 1 (121) 
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! (B + F)? y e x + P 2 + (qDP 2 * i F )*y 9 y = P 2 t 12 ^ 

Equations (121) and (124) yield 

( 125 ) 

In connection with equation ( 100 ) , equations ( 85 ), (86), and ( 87 ) ^PP^y 
here as veil, so that equation (2) becomes, analogous to equation (ott), 

P cr = P 0 + ^ ( 0 3cPUc + V^) 



(126) 
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where Pq can usually "be neglected. For simply supported plates, from, 
references 11 and 12, 


p 0 - 


2x%I f 


+ 2(B 
P 2 


+ 2F) + Dp" 


where (fig. 3) 

I f = h 5 /l2 


(127) 


( 128 ) 


For clamped plates the deflections w^y and Wgy of the Y-strips are 
shown in figure 15 . Hence only for the deflection wjy do the faces 
act as clamped, while with the deflection Wgy they practically act as 

simply supported so that in accordance with item (5) in the section 
"Description of Method" 


p o “ p oi + p 02 ~ e y p oi + ( 1 " e y) p 02 ( 12 9 ) 

where Pq^ is identical to Pq from equation ( 127 ) and, since the clamped 
edges p = 2.5 and q = 5 (refs. 11 and 12), 


2jc 2 E, 


01 


*f 


A + 2.5(B + 

P 2 


2F) + 5DP C 


(150) 


Equation (129) is slightly conservative, since for case 2 the break in 
the faces at y = 0 and b requires extra energy. Hence, for simply 
supported unloaded edges, from equation ( 126 ) using equations ( 113 ) and 
(127) with p = 2 and q = 1, 


, sA B I f 
or " b £ 


A + 2(B + 2F) + Dp 2 


P 




ffil.Cy) + (B + 2F) 


e x + [(B + 2 F) + Dp^l 8 y l (131) 
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For clamped sandwich plates, from equations (113), (126), (127), (129), 
and (130), with, p = 2.5 and q = 5> 



(132) 

where 0 X and 0y are given by equation ( 125) • Besides from the 
formulas in references 33 and 34 A, B, D, and F can he read directly 
from charts in reference 39 as functions of the secant and tangent moduli 
for the actual "buckling stress c cr = P cr /2h. 

In the elastic range, from references 33 or 34, 

A = B + 2F * D = l/(l - V 2 ) (133) 

so that for a simply supported plate equation (131) reduces to 



(134) 

where 0 X and 0 V from equation (125) with p = 2 and q = 1 reduce 
to 
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( 135 ) 


If the proper flexural rigidity of the faces is neglected (l, = 0 ) and 
P cr is written in the form ' ' 

it 2 E„I„ 

P=r * k s —4"^ (136) 

h^ 

k s from equation (136) can he shown to he identical to k s from refer- 
ence 32. With p = 0, since in equation (135) ? x = (l - v2)r x n.nrl 
?y = (l - v2 ) r y> these equations for 8 X , Qy, and P cr reduce to 
equations (83) and ( 89 ) j so that also these equations are exact. 

For a clamped sandwich plate in the elastic range equation (132) 
reduces to 



31 E sls 1 + 

(l - v 2 )h 2 l _ 


( 137 ) 
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where, with p = 2.5 and q. = 5 > from equation (125) 

i‘ . 5 » 8 * - v) f, - 0.625 r^, I, 


H5 


l 


1 * (q t ihi + 0.625(1 - v) 


» a *[ 5e ‘ * °-® !i - v i]r^) -sfrSj 


1 - V' 


p 2 + 


1 + 


1 — y2 ^ + 0.625(1 - v) 

p 2 


[ ~ 1 + i ( l- 2 v2 ) f 1 + 0 , 6 ^_ d ; - 0.625 ji - 




Ip 2 + 


o "I r_ 

5P + 0.625(1 - v) — £ 
Ji - 


6.25?. 


16(1 - V 2 ) 


l (158) 


For isotropic core, r x = r y = (l - v 2 )r where r is given by equa- 
tion (78) or (8l) with G x = Gy = G c or 


r = 


* 2 e s i s 


(l - v 2 )b 2 G c (t + h)‘ 


(139) 


Here G c is the modulus of rigidity of the core. Further H = 0. Thus 
equations (137) and (1-38) reduce to 


P 


cr 




(1 - v^)l 


P + 2 + \ V + 




+ 


11 

(l - v 2 )b 2 


(p + 1 - 2 5) 0 * + 


(1.25 + 5f3 2 )e y 


(i4o) 
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where 


p 2 + 

>75P 2 + 0.625(1 - v)(l + p 2 ) 

u ] 

P 2 + 

1 + 5P 2 + 0 . 625(1 - v 

)(l+ p 2 )jr + 

3A375 + 0.625(1 - v)^-L + 2.5 + 5P 2 j 

r 2 


> (aJft) 


P 2 - 

[_Q.25 - 0.625(1 - V)(l + 0 2 )_ 

. r 

P 2 + 

1 + 5P 2 + 0.625(1 - v) (1 + p 2 ) 

r + 

3 A 375 + 0 . 625(1 - v)| 

^ * ^-5 ♦ ® 2 ) 

*2 


It seems that the maximum discrepancy between the simple formulas 
for isotropic- core sandwich plates derived in reference 11 and the exact 
calculation occurs for the present case of clamped edges. From refer- 
ence 28 both lead to the same buckling loads for r = 0 and r ^ 1 and 
the maximum discrepancy of 7 percent occurs if r is abcrut 0 . 25 . From 
table 1 of reference 28 where, as usually can be done, the proper rigidi- 
ties of the faces have been neglected, a minimum value of the buckling 
stress coefficient 


(lte) 

n^N a rt“E s I g 

occurs with r = 0.25 for P = 0.6 where k = 2.88, while from refer- 
ence 11 for the same case k = 5*08. Also equations (l4o) and (l4l) 
with If = 0 for r = 0 and r ^ 1 lead to the same values of k as 
these given in reference 28. For the case of greatest discrepancy 
between references 11 and 28, with r = 0.25, P = 0.6, and v = 0.333., 
they yield 0 X = 0.524, 0y = 0.274, and k = 2.94, or only 2 percent 
more than the exact k value from reference 28. On the other hand, 
‘though more accurate, equations (l4o) and (l4l) are more intricate than 
the equations derived in reference 11 . Their advantage with respect to 

those in reference 28 is that an explicit formula is obtained for P 
. cr 

or k. 

* - 

For sandwich plates with orthotropic core, where as stated in refer- 
ence l 4 „the method of references 11 and 12 cannot be used if the 
ratio Gy/G x differs too much from 1, the present formulas lead to very 

accurate' results. It was shown in the foregoing discussion that for 
simply supported plates they are exact. For clamped plates, in the 
limiting case where r x = 0 and with p. = 0 and If = 0, equations (137) 

and (138) reduce to 
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P 


cr 


jr%! I 
B E 

(1 - v 2 )b 2 



+ 



0 X 


+ 



(1^3) 


-v 

e x = 1 


f. 2 - 

[o. 625 r y /(l . 

- v) 



r 

r 5 ^ + 

L_ 

0.625(1 - v)j 

IV 

t 1 - v2 )) 


omj 


From these formulas, for example, with 0 = 1, r y = 0.25, 1> and <», 
and v = 0 . 33 , one finds that 


*8 = 



* 2 e s i s 


( 1 ^ 5 ) 


is equal to 4.67, 2.59, and 1.45, respectively. From figure 4(a) of 
reference 52 for 0=1 and the same values of r y , k s is 4.55, 2.51* 
and 1.4o, respectively, so that the discrepancies are not more than 
3 percent. 

This shows that also for sandwich plates with orthotropic core the 
method of split rigidities leads to accurate results. If necessary, 
similar formulas can he derived for other boundary conditions. 


Buckling of Homogeneous Plates Under Nonhomogeneous 
Stress Distribution 

Item (6) in the section "Description of Method" can be used to calcu- 
late deflections of columns and plates with initial crookedness, as was 
done in footnote 2 of reference 1 and in reference 15, respectively. It 
can also be used to find the axial resistance of initially flat plates 
under nonhomogeneous stress distribution, as was done in references 15 
and 17 for plates in the postbuckling stage.. The same method can be 
applied for calculating the critical load for incipient buckling of plates 
under initially nonhomogeneous stresses. In all cases the restraining 
or internal action of the plate is derived from the case of uniformly 
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distributed stress. Hence the condition for the accuracy of this method 
is that the shape of the deflection surface for the nonhomogeneous stress 
distribution considered is practically similar to that under unif orm 
stress. 

Such a case occurs, for example, with buckling of a long simply 
supported flange under linearly distributed stresses 0 X (fig. l6). In 
a similar way as under uniform stresses a x the bending moments My 
(and M x ) will be relatively very small so that the plate can be assumed 
to remain straight in the lateral direction or 


w = Ky sin | x (l46) 

For uniformly distributed compressive stress with Poisson's ratio equal 
to 0.3, from reference 40 


"cr - * ^ 

The external and internal actions that are compared (item ( 2 )) are the 
moments of the deflecting and restraining forces with respect to the 
hinged edge. For uniform, compression the deflecting forces acting upon 

a small element t dx dy are -tcr cr (d^w/dx 2 ) dx dy, so that from 
equation (l46) they can be expressed as 


b 2 !: 


(iVf) 


D = C 1 cr cr w dx dy = Ca cr y dy (l48) 

where C is independent of y. Hence the external moment exerted by 
the forces D acting upon a cross strip b dx about the hinged edge 
y = 0 is 


r r u * 

Me U = / Dy dy = 00^ / y 2 dy = J Ca cr b^ (l^9) 

J 0 Jo 5 

Then, from item (6), the moment exerted about the hinged edge by the 
restraining forces acting upon the strip must be 


M i = Meu " k Ccr cr t ' 


(150) 
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Since under nonuniform compressive stresses (fig. l6) equation ( 1^+6 ) 
remains valid, also in that case is given hy equation (150). With 

compressive stresses 


°x = 



(151) 


the deflecting forces are, analogous to equation (l^8), D = Ca x y dy, 
so that they exert a moment about the hinged edge of 


Me = F Dy dy = 

Jo 




dy = Ca e 




dy 



Ca e b 5 


(152) 


Since M e should be equal to M^, from equations ( 1^7 ) , ( 150 ) , and (152) 
the critical maximum edge stress a Q is 



1-70 

h ~ T b^ 


(155) 


For t = 1 and 2 this yields k values of O .568 and O. 85 O, respectively, 
in accordance with values given in reference 40. If the maximum stress a g 
occurs at the hinged edge (fig. 17 ), so that 


= (?p1- 



(15*0 


and t varies between 0 and 1, the external moment is 



( 155 ) 
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s o that • from equations ( 1^7 ) > ( 150 ) > and- ( 155 ) 



J cr 


= 1»70 
h - 3t h 2 t 


(156) 


For r = 1 one obtains k = 1*70, in accordance with reference 40. 


Ultimate Load of Plates Under Compression 

In reference 15 item (6) (section "Description of Method.”) was used 
for determining the postbuckling behavior of a simply supported plate, 
the unloaded edges of which are held straight in the plane of the plate, 
assuming the plate to remain elastic. This method will here be extended 
to the plastic range in order to find the ultimate load of such a plate 
(ref. 4l) . 

It was shown in reference 15 that up to deflections of the order 
of the plate thickness t practical exact results can be obtained by 
assuming the shape of the deflection surface of the plate to remain 
similar to that at incipient buckling, that is, 

w = w m sin £ x sin £ y ( 157 ) 

where a = b. Considering the deflection to be developed in a Fourier 
series, the distribution of the total direct stress cr x in the post- 

buckling range, as sketched in figure 18 , will tend to superimpose on the 

deflection from equation ( 157 ) partial deflections Aw = w-^ sin ^ x sin ^ y, 

for which the individual equally distributed buckling stress a x is 

25 times that for the mode of equation ( 157 ) • The stresses ay (fig. 18) 

superimpose deflections w^ sin x sin ^ y and w^ sin 22 x sin ^ y, 

with individual equally distributed buckling stresses <y that are 25 

and 9 times, respectively, those connected with equation ( 157 ) • Hence 
the required deflecting forces for bending in these modes are much higher 
than for bending in the mode of equation ( 157 ) > which makes the contri- 
bution of these modes relatively small. In the elastic range the mem- 
brane stresses a ^ and a in figure 18 are (ref. 15) 
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(158) 


In reference 4l it was shown that it is sufficiently accurate to assume 
the same distribution in the plastic range. Using a method that differs 
somewhat from that used in reference 15 the uniform stress distribution 
( ff x)u "that is equivalent to the membrane stresses cr^ can be calcu- 
lated as follows: If a x = (cr x )u does not vary with y the work done 

by these compressive stresses per buckle with deflections w is from 
equation (157) 


°Tny - <%xm ^ y 

°my = a mym sin2 l x 


Vxa ^/ 0 f Q W u (^) ^ ^ 2 W u a2 ** & 4 (159) 

On the other hand the tensile stresses from equation (158) would 

exert a work 



(160) 

From equations (159) and (160) the membrane stresses 0^ are equivalent 
to equally distributed compressive stresses 



_ 5 


u 


mxm 


The work done by uniform compressive stresses f< 3 y \ u I s 


(161) 


'yu 


■t/o7o 



^ (162) 
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so that from equations (159) and (162) an equally distrihuted stress cr y 
is equivalent to a stress y 

( CT x)u = ^ ( ff y)u = P (°y)u (^ 3 ) 

b 

Hence, from equations (l 6 l) and ( 163) the total load on the plate as 
shown in figure 18 is equivalent to a uniform load 

(°x)u = a xe ” ^ °inxm + P jg CT iny m ” 4 a mymJ (l 64 ) 

The membrane stresses are caused by the increase in length of t he 

middle strips GH due to the deflection w which amounts to 


6 



( 165 ) 


so that the difference between the average direct strains in the strips AD 
and GH is 


^ - t (5 s ) 2 <«> 

Assuming, as proved to be true in the elastic range (ref. 15 ), that c r x 

does not vary with x and that Poisson's ratio does not influence the 
stress distribution so that it can be assumed to be zero, this gives a 
stress 


a = E e - rt2 e r m | 
°mxm Ji Bx e xm ~ ^ ^sxl~J 


(167) 


where E sx is the secant modulus to be applied. Similarly 


mym 


* 2 _ N 1 


' J sy (L ym IT E sy^ b 


(168) 
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bo that 


^mym P eI? a mxm ( 169 ) 

sx 

At incipient bu ckling the restraining forces offered by the bending 
rigidity of the plate are just in equilibrium with the deflecting forces 
caused by the compressive stresses cr x = T)a cr , where cr^ is the elastic 

buckling stress of a simply supported plate (ref. 2 6, p. 529 ) 


a cr - (f» * I)' ^ ^ (1T °) 

and T] is the plastic reduction factor. The deflecting forces are pro- 
portional to the center deflection w m of the plate. Since the restraining 
forces are due only to the bending stresses in the plate (the membrane 
forces are included in the loadings <Jx s.nd cfy) in the elastic range 

they are also proportional to the deflection w m . Therefore, with finite 
deflections the plate will be able to resist stresses a x and cr y 
(including the membrane stresses) that are equivalent to uniform stresses 

(°x)u = T l ff cr (iTl) 

where rj is the plastic reduction factor of the plate flexural rigidity 
at the actual finite deflection w m . Since the actual loading of the 

plate is equivalent to (cr^ u from equation (16*0, combination of equa- 
tions (l 64 ), (169), and (171) gives 


1 a 1 * „ 


a xe “ % a mxm 4 ^ CTmxm 

SX 


(172) 


or 


H g xe ~ ^ P cr) 
' W " 3 ♦ (pV/Esx) 


(175) 
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( 158 ) the average stress, that is, the 


° rny <3y - ff xe " 2 (174) 

Apparently the postbuckling stress depends on (3 = a/h. In general, 

the value of J3 corresponding to a minimum value of cf^p will decrease 
as the deflection v m increases. However, as explained in reference 4l, 

the plate will not he free to change its wave length continuously because 
at a certain wave length the buckles will "freeze." This process can be 
illustrated by considering a very long plate for which any wave length 
would satisfy the boundary conditions . For such a plate the plastic 
deformation due to compressive and bending stresses will cause permanent 
buckles so that after a certain stress has been reached the wave length 
becomes fixed. From available test results it appears that the wave 
length that establishes itself at incipient buckling is maintained in the 
postbuckling range. The only plates that exhibit a postbuckling stress 
are those that buckle in the elastic range, where at incipient buckling 
P = a/b = 1. Hence in equation (173) ^ may be assumed to be 1. In 

that case E sy ^E sx may also be assumed to be 1 . 

From references 25 and 42 for simply supported plates with edges 
that are not held straight the ultimate load is a function of the yield 
stress ffy S and the critical stress cr^ . The (average) ult im ate stress 
was found experimentally as 

°ult = ( l 3 e/ t ) a ys = \j 0 cr°ys “ 0 . 25 cy cr (175) 

Ass uming that in the present case, where the edges are held straight, 
the ultimate load is reached when the edge stresses cr xe reach the yield 

stress, from equation (173) with P = E S y/E sx = 1 and b xe = 0y S , 

°mxm = °ys ~ 4°cr (176) 

so that the ultimate stress is, from equation (174) with and 

0 xe = °ys* 


Hence from figure 18 and equation 
postbuckling stress 0 ^, is 



a ult “ a ys 


1 a 

2 mxm 


(177) 
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In reference 4l it was shown how for a given case tj can he calculated 
from the plastic deformations of the plate. Since this calculation is 
very involved a direct formula was derived for tj, which for the case 
considered was shown to he in excellent agreement with the more exact 
calculation. From reference 17 for edges that are not held straight in 
the plane of the plate, to which equation ( 175 ) refers, the membrane 
stresses c^ and 0 ^ are equivalent to a uniform stress 

( G x)u = -O-^mxm UT8) 

Also in that case the equivalent loading of the plate should he equal 
to Vcr* With- edge stresses cTy g this gives 

0y s - 0 . 848o - mxm = TjCgj, ( 179 ) 

or 

°mxm = l*l 80 ys - 1. l8qcr cr ( 180 ) 

Equation ( 17*0 also applies here so that, with = a ul -(- and cr xe =» 0 ys , 

using equation ( 180 ) gives 


°ult “ a ys “ 2 a nixm = ^ • ^°ys + 0.56*i G cr (l8l) 

and 

n = 1.7 5 ^ - 0.7 5 s ( l82 ) 

°cr °cr 

Insertion of from equation (175) yields 

1 

i = i.7&^) - 0A25 - 0.7 ^ (185) 

\ u cr/ cr 

Since 1 } thus is based on test results for relatively large deflections, 

it automatically also contains a correction on the use of equation (157) 

in that range. In order to apply this result to the present case with 

edges that are held straight, it should first he determined on which 

variables n depends . It Is dependent on the direct stresses a v and a 

a y 
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that vary with (fig. 18) and on the "bending strains measured "by 

the strain difference "between "both faces of the plate. At the center of 
the "buckle this difference is 


Ac _ j- _ jt£ „ + 

•xb ~ ^2 a 2 


On the other hand, from equation (l 66), 


v - ‘ f (%») 1/2 ' 


"m 


2a/ g mxm ^ 
* l^sx 


1/2 


so that from equation (184), with a = b, 


A 2irfc /°mxm) 
Ae *b = — i: 


1/2 


J sx 


From equation (170), with 3 = a/b = 1, 
so that equation (l86) becomes 


± = 0.525 fej 


\l/2 

Ae x -b = 3.3 


' g cr g mxml 

e e sx ; 


(184) 


(185) 


(186) 


(187) 


(188) 


Hence, since t] depends on 0^ and Ae^, it obviously can be con- 
sidered as a function of and cr cr . From equation (179) 

is = 0.848 is + , (189) 

°cr °cr 

Inserting this into equation (185) gives T) as a function of 

CTqj. 
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n = 0.25 - 0.35 


J mvm 


'cr 


+ 0.71| 


' g mxm \ 
?cr ) 


1/2 


(190) 


With this value of 7} for the present case, from equation (176), 

fa^ — N 1 / 2 


cr 


S = 0.25 + 0.35 — - 0.71 P™ 1 


cr 


'cr 


'cr 


(191) 


from which 


mxm _ 
^cr 


-0.5^5 + fl-5^ ^ - O.O865W 2 


(192) 


After calculating 0,3™ from equation ( 192 ) , follows from equa- 

tion (177). 

In figure 19 a u it/°cr froin equation (175) has "been plotted as a 
function of d C r/ a ys a8 the dashed curve (l) . In the same figure perti- 
nent test results for simply supported plates with edges that are not 
held straight in their plane are indicated by squares, as plotted from 
reference 43 • The squares refer to square-tube tests by Needham (ref. 44) 
and the diamonds refer to V-groove tests by Anderson and Anderson 
(ref. 45) which are also shown in reference 43. This shows that formula 
(175), determined by Winter from tests with light-gage steel (ref. 42) 
involving plates that are not held straight in their plane, leads to 
excellent results. Hence, the formula for rj, derived from it and pre- 
sented herein as equation (190), can be considered as well based. There- 
fore the combination of equations (177) and (192) can be expected to 
yield satisfactory results for plates with simply supported edges that 
are held straight in the plane of the plate^ In figure 19 a uit/ a cr 
from these equations has been plotted against a C r /°ys aB "the solid 
curve (2); the test results from Botman (ref. 46), as given in refer- 
ence 43, are indicated by triangles. These tests were done with three- 
bay plates between knife edges so that the edges were not held completely 
straight. The test results should therefore be located between the 
curves (T) and (g) , which indeed occurs. The theoretical results from 
reference 43 Tor plates with edges held straight are indicated by circles. 
These are a little above curve (2) . This may be because in reference "43 
Poisson's ratio was assumed to be 0.5 for plastic as well as elastic 
deformations which increase the plate rigidity above its actual value. 
Hence equations (177) and. (192) yield reliable results. 



NACA TN 4085 


61 


Formulas based on the same principles have been derived in refer- 
ence 4-1 for the ultimate load of simply supported and clamped plates 
under the additional influence of thermal stresses. Also formulas were 
derived for the equally distributed stresses ( ff x)u ( a y) u '^ La ^ i 

are equivalent to arbitrarily distributed thermal stresses c x and cy 
as far as critical as well as ultimate loads are concerned. 


PROBLEMS WHERE COMPONENT CASES HAVE DIFFERENT BOUNDARY CONDITIONS 


Problems where the component cases have different boundary condi- 
tions have been dealt with in references 13, l6, 17, 18, and 19 . As 
pointed out under item. (6) these problems lead to more intricate relations 
between P cr and P-j_, Pg, . . . P n than those where the boundary condi- 
tions for all cases are the same. A simple example is' the problem of 
deter mining the bu ckling stress of an I-section about its minor axis. 

The buckling deflection of the cross section is shown in figure 20(a) and 
is split into case 1, a deflection as a column without distortion of the 
cross section (fig. 20(b)) and case 2, a distortion of the cross section 
(fig. 20(c)). This problem was dealt with extensively in reference 13 
where the deflecting and restraining forces were chosen as the external 
and internal actions. This method was sufficiently accurate because the 
decrease of the column buckling stress a - (case l) or plate buckling 
stress cfg (case 2) from the interaction of cases 1 and 2 was only a 
small fraction of the actual buckling stress. 

Somewhat more accuracy can be obtained by comparing the work done 
by deflecting and that done by restraining forces. This method was used 
in references 17, 18, and 19 . It will be demonstrated here in deriving 
the buckling stress of stringer panels. 

The stringer panels considered are supported by stiffeners of equal 
cross section and located at equal distances b (fig. 21). The half 
wave length of buckling in the X-direction is denoted as a (fig. 22(b)). 
If no stiffeners were present, the buckling deflection of the plate would 
vary sinusoidally in the X- and Y-directions . If stiffeners are present 
and the deflection w-^ (fig. 22(a)) which determines the deflection of 

the stiffeners varies in the same way, the additional deflecting forces 
-ta x (8^w/8x 2 ) per unit plate surface, caused by the additional compres- 
sive stresses c x which the panel can sustain due to the stiffeners, 
will as a whole vary sinusoidally too. Hence, equilibrium is possible 
if the restraining forces El(d\r/dxO per unit length supplied by the 
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stiffeners also vary sinusoidally in both directions. This is indeed 
the case if w^ varies sinusoidally. Hence, it may obviously be assumed 
that 


W 1 = v lm sin ? Sin S = Y 1 sin ^ 


(195) 


where 

\ = it/a (194-) 

In a similar way, as was done in reference 13 (fig. 20 ) for deter- 
mining the interaction between column and plate buckling, the buckling 
deflection of the panel (fig. 22(a)) could now be split into two caBes, 
and the buckling stress 0 cr expressed in terms of the individual 
buckling stresses for these cases. In the present case a slightly 
different approach is used by first splitting off the buckling stress a Q 
of the unstiffened sheet, caused by its buckling into the shape w^ from 

equation (193) and denoted as case 0, as indicated by the first method 
under item (5)« From page 329 of reference 26, 


where 



_a_ nb\ 2 l\ 2 

nb a] n 2 b 2 t “ ( n 2 p) ^ 


^ = a/b 


(195) 


(196) 


and N is the flexural rigidity of the sheet. The resulting buckling 
deflection of the reduced structure with buckling stress a r = cr cr - Oq, 

where consequently the partial buckling stress cTq caused by the rigidity 
of the sheet against bending into the shape w^_ is neglected, is then 
split into two individual cases. Case 1 is the deflection W]_, shown 
separately in figure 22(c). Case 2 is the extra deflection w 2 , shown 

separately in figure 22(d), which is caused by the deflecting forces acting 
on the sheet and created by the partial buckling stresses a r = a cr - oq 

that cannot be taken by the unstiffened sheet. These deflecting forces 
are consequently resisted by the stiffeners, which cause the bending of 
the sheet into shape w 2 like a laterally loaded continuous plate on 
several supports. 
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Hence, the actual buckling stress of the panel is 

^cr = ff 0 + (197) 

•where a r has to be expressed in the buckling stresses and a 2 for 

the individual cases 1 and 2, respectively. 

Since the flexural rigidity of the sheet has to be neglected here, 
in case 1 the buc kli ng stress is due "to the flexural rigidity of the 

stiffeners alone. The critical thrust of the single stiffeners, with 
effective moment of inertia I, is 

P x = ir 2 El/a 2 = A 2 EI (198) 

Hence, the restraining force acting on an element of length dx of 

a stiffener if bent into a half sine wave is equal to the deflecting 
force = -P^ (d 2 w-j_/ dx 2 ^ dx or, from equation ( 193 ) ) 

R 1 = D 1 = ^ 2p l w i (199) 

At incipient buckling of the complete panel the compression ^ B ±, a cx dn 
the stiffener itself, with cross section A s ^, causes a deflecting force 

D s = -A s tcr cr ^wx/Bx 2 ) dx = X 2 A B tff cr wi dx (200) 

so that the restraining force exerted by a stiffener on the sheet is only 

V = R 1 - D st = A 2 (P x - A Bt cr cr )w 1 dx (201) 

Therefore, by comparing equations (199) s-nd (20l), the restraining force 
exerted by a stiffener on the sheet is equal to that exerted by a stiffener 
with a zero cross section and a critical thrust 

p l’ = P 1 - A st*cr “ ^ - Ast^cr (202) 

The buckling stress for buckling in a shape like that in figure 22(d) 
may be written in general as (ref. 25) 
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’2-^ = (f +p + #2 )^ (£05> 

where the computation for p and q is given in the appendix. However, 
since only the reduced case is considered, where the part Cq from equa- 
tion (195) for buckling in a shape like w^_ must he ignored, the critical 
stress for the present case 2 is 






b^t 


(204) 


In order to express the reduced critical stress a T of the combined 
structure in terms of P-j_' and bg' from equations (202) and (204), the 

work done by the deflecting and restraining forces will be compared. The 
deflections w^ and W2 in figure 22 both vary sinusoidally in the 

X-direction (fig. 22(h)) so that, similar to equation (193), 


v 2 = Y 2 Bin Ax 


(205) 


where Y 2 is a function of y alone. Hence, at the reduced buckling 

stress the deflecting force in the Z-direction (fig. 22(a)) exerted on 
a sheet element t dx dy as a consequence of the deflection w = v-^ + Vg 

is -tcr r (d 2 w/dx 2 ) dx dy or, from equations (193) a-^d (205), 

D = A ta^w^ + w 2 ^ dx dy (206) 

Since from figure 22(a) the deflection w 2 does not involve any 
deflection of the stiffeners, it is evident that the work done by the 
deflecting forces for the deflection Wg is taken entirely by bending 

of the sheet in the shape w 2 * Furthermore, since the reduced case is 

considered, where the flexural rigidity of the sheet is ignored for the 
shape w 1 , the work done by the deflecting forces for the deflection w-j_ 

is taken by bending of the stiffeners alone. 

The deflection w^ will be considered first. Denoting the deflec- 
tions of the stiffeners as 


w st - w stm 


sin Ax 


(207) 
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where w s ^ = W]_, from equations (201) and (202) the restraining force 
exerted by a stiffener per length dx is 

%' = ^lXt to ( 208 ) 

From the foregoing discussion the work done by the deflecting forces D 
from equation ( 206 ) for the deflection w-j_ has to be equal to the total 
negative work done by the restraining forces R-^ 1 from equation ( 208 ) 
furnished by the stiffeners. Hence, 


zJT 1 ** 1 ^ ^ = ijj R i v i ^ dy 

or, using equations ( 206 ) and ( 208 ) and dropping the terms A^/2, 

tOp JJ' Vj 2 dx dy + tcr r JJ ' v^Wg dx dy = ^ w st 2 to ( 209 ) 


where the double integrals have to be extended over the entire sheet with 
dimensions a times hb and the summation sign refers to the (n - 1) 
stiffeners. Equation ( 209 ) may be written in the form 

nbtcr r (w-^ + cpwgm) = (210) 

where w-] m and w^ are indicated in figure 22(a) and where 


JJ .jVg dx dy 


© = Js 

w 2m rr T.r ^ 


( 211 ) 


JJ w 1 dx dy 

and, since w^ = w 1 (fig. 22(a)), using equation (195) gives 


i=n-l 


t = 


^ l L'* ^ ^l^t 2 


nb 


Z sin2 ^ 


ib 


i=l 


^dxdy J o v x 2 dy J* 


= n 


( 212 ) 


. 2 « , 
sin — y dy 
nb 
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Hence equation (210) transforms to 


(a r - a^w^ + q><V*2ra = 0 

where, from equation (202), 



kct, 


cr 


(215) 


(214) 


rr - V USS _ ^1 _ At = 1 El A 

1 " 1 b 2t bt " a 2bt i 5 m b 2 ! 


(215) 


K 



( 216 ) 


Replacing cr-^' in equation (213) hy the last term of equation ( 2l4 ) , in 

which again 0 cr is expressed according to equation ( 197 ) > equation ( 213 ) 
becomes 


(1 + K)ff r - - K0q) 


v Im + c P 0 r w 2m = 0 


(217) 


In order to obtain a second relation between w lTn and Wp^ the work 
done for the deflection w 2 is considered. In case 2 buckling occurs 
at the stress Og' from equation ( 204 ) , so that the restraining force 
R 2 acting on an element t dx dy is equal to the deflecting 
force D £ = -tcr 2 ' ^8 2 w 2 /dx 2 ) dx dy = A 2 ta 2 ’w 2 dx dy. This restraining 
force R 2 depends on the deflection w 2 alone and not on the compressive 

stress so that, with the assumption that in the combined case the deflec- 
tion w 2 has the same shape as in case 2, also in the combined case the 

restraining force is 


R 2 = A tffn'Vp dx dy 


(218) 
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Equating the work done for the deflection Wg "by the deflecting forces D 
from equation ( 2 06 ) to the negative work done by the restraining forces R 2 
from equation (2l8) gives 


iJJ I>V2 ^ ^ = ijf b 2 w 2 ^ ^ 

Using equations (206) and (218) and dropping the terms A^/ 2 , this 
becomes 

dx dy + ta r JJ' v 2~ ^ ^ = ta 2' Jf 4x1 dy C 21 ^) 


or 


where 


7^r w lm + ( CT r “ 0 2') w 2m = 0 C 220 ) 


7 = 


w, 


2m 


w. 


Im 


// 

// 


w^wg dx dy 


2 

Wg dx dy 


( 221 ) 


The buckling condition is obtained from the two homogeneous equa- 
tions ( 217 ) and (220) by equating the denominator determinant to zero, 
which yields a T , so that from equation ( 197 ) 

0 cr = a 0 + CT r 


cr-^ - kctq + (l+ / 

r > 

Hi - 79 + k) (a x - K0 O ) eg’ 

iM 

ff o + 2(1-7*?+ k ) r 1 

[ 0 1- Kff 0 + (l+ Kjag'j 2 
^ j 

/ 


( 222 ) 
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or 


„ . k 4 1 

° r t 2 t 


(223) 


where 


k cr - k 0 + 2(1 + k) ( k l ‘ Kk o) 1 1 " \ 


^ _ 4(1 - yq? + *)ft 
1 + n(i + k )] 2 


l/2\ 


ft = ^2 ' j (-^i " Kjs o) 

and from equations ( 195 ) # ( 204 ) , (215)# and (216) 

*-(***)’ 


in ,111 

*1 „2 hN 


l\ q 2 


V " p - -f + «'iP 


A. 


K = 


St 

ht 


(224) 

(225) 


(226) 


J 


The values of p, q# and 79 for the cases of 1, 2# 3# and an infinite 
number of stiffeners are determined in the appendix and are assembled 
in table 2. If Is^ = 0 and k = 0 so that no stiffeners are present# 


equation (224) reduces to k^ = kg# as it should be. 


thus ft 
becomes 


are infinite, so that no interaction occurs# 


if ( is 

ln2 3j p 2 bN 


If k2 1 and 
equation (224) 


( 227 ) 
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in accor dan ce with, equation (l) of reference 47* Equation (227) niay 
tie used if {3 > 8. The largest interaction occurs for small values of (3. 
In table 3 the results from equation ( 224) for {3 = 1 are compared with 
those scaled from the graphs of reference 47- The values from refer- 
ence 47 for k cr > 4 refer to the extension of the El/(bN) curves. 

The only important discrepancies occur for n = 3* El/(bN) = 0, and 
k = 0 , where the value from reference 47 is placed within parentheses, 
since here indeed from equations (224) and ( 226 ) k cr = Eq = 100/81 = 1.235 

so that the graph seems to be inaccurate at that point. 

The interaction is the smaller the larger the value of ft is from 
equation (225). In a similar way as was done with equation ( 69 ) of refer- 
ence 13 , for not-too-small values of 2 equation (224) may be simplified 
by observing that the term in the large parentheses may be written as 


1 - (1 - e ) 1 / 2 = l- (l-^-^ 2 )=|e(l + ie) ( 228 ) 

since in that case e is small with respect to unity so that equa- 
tion (224) transforms to 




k cr k 0 + 1 + n(l + k) 


1 + (1 ~ 79 + *)» 

[l + 0(1 + k)] : 


'(kq - Kk o) (229) 


For still larger values of ft, the term 1 in the denominator of the 
fraction within the large brackets may be neglected, so that equation ( 229 ) 
becomes 


cr 


~ *0 


ft + 


1 - 79 + k 
(1 + k ) 2 


- >ckQ 

1 + ft(l + k) 


( 230 ) 


Equations ( 229 ) and (230) are sufficiently accurate if ft is larger 
than 1 and 4, respectively. 


Values of k Qr may be calculated directly from equations (224), 

( 229 ) , or ( 230 ) . They have to be calculated for half wave lengths 
a = if m, where l is the simply supported length of the panel in the 
X-direction (fig. l) and m is an integer (for only an infini te number 
of stiffeners can a be assumed to be equal to l). The lowest value 
of k cr so obtained has to be inserted into equation (223). The 

effective moment of inertia I of the stiffeners to be used in 
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equations (226) for calculating k^ has to be determined as shown in 

reference 20, where also a formula for the buckling stress coefficient 
for forced crippling Is derived from equation (224). 

This example shows that the method Is applicable and gives accurate 
results for very intricate problems as is shown also by its application 
in references 13 and 17. The present problem was also worked out for 
stresses a acting simultaneously in a direction perpendicular to the 

stringers as may be caused by thermal stresses, but it was thought that 
omission of the stresses Cy gives a clearer picture of the method. 

It should be noted that in reference 13, where the deflecting and 
restraining forces acting upon an arbitrary small element were chosen 
as the external and internal actions, these forces were expressed in 
terms of the maximum deflections w, w^, and W£ to which they were 
proportional. In the present case, and also in reference 17, where the 
work done by the deflecting and restraining forces was considered, the 
forces acting upon all elements must be taken into account. These forces 
cannot be expressed in terms of the maximum deflections, since for dif- 
ferent elements different proportionality factors apply. Therefore they 
were expressed in terms of the local deflections w, w-^, and Wg of 

the individual elements, which expression in all cases considered was 
facilitated by the sinusoidal variation of the deflections in the X- 
direction by which the deflecting forces -tcr x d^/dx 2 are proportional 
to w. 

Since references 13 and 17, where the method was applied to several 
other cases with different boundary conditions of the component modes, 
are directly available, It is thought unnecessary to give other examples. 
The reader therefore is requested to consider references 13 and 17 as 
part of this section. 


Cornell University, 

Ithaca, N. Y., June 28, 


1956. 
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APPENDIX 


COMPUTATION OF p, q, AND 79 

Computation of p and q in Equation (203) 

For n = 2 .- For plates with one stiffener (fig. 23(a)) in case 2 
(fig. 22 (d)) each plate of width b can be considered as simply supported 
at one unloaded side and clamped at the other one. From reference 48 
for the clamped simply supported case p = 2.27 and q = 2.45. 

For n = 3 «- For plates with two stiffeners (fig. 23(b)) plate AB 
is simply supported at edge A while at B it is subjected to a moment 
that varies sinusoidally in the X-direction. Plate BC is subjected to 
equal and opposite sinusoidal moments at edges B and C. The buckling 
coefficient kg can be determined by requiring that the spring constant 

against rotation of plate AB along B is equal and opposite to that of BC 
along B. Spring constants for these cases or rather the specific angle 
rotations were derived in reference 49, where they are given by equa- 
tions (48), (52), (53)? and (54). The same results were obtained in 
reference 50 and are presented in tables in reference 51 , where the 

stiffnesses S 11 and S 1 ^ are the quarter spring constants referring 
to the present plates AB and BC in B. Hence, k^ is determined by the 

condition S ’*- 1 = -S^. It appears that the minimum value of k^ occurs 

for p = a/b = O .72 and is equal to 6 . 15 . From equation (6.4.1) of 
reference 25 this minimum is 

Wmin = P + 2 ^ (Al) 

Inserting the value 6.15 into equation (Al) and kg = 6.935.? obtained for 
P = 1, into equation (203), one obtains p = 2.36 and q = 3 >58. 

For n = 4 . - In the case of n = 4, (fig. 23(c)), from symmetry 
plate BC can be considered as clamped at edge C so that, with the nota- 
tions of reference 51 , the condition at B is S-^ = -S, yielding a 
minimum value of kg = 6.47 and, for p = 0*9? kg = 7* 00. From equa- 
tions (Al) and (203) this results in values of p = 2.37 and q = 4.20. 

For n = 00 . - For the case of n = <» all plates can be assumed to 
be clamped at both edges (fig. 23 (d)) so that, from reference 25 , p = 2.5 
and q = 5- Values of p and q are assembled in table 2. 
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Confutation of 

From equations (21l), (193).> and- (205)> 

// v-,Wo dx dy / Y-,Yp dy 

_ JJ 12 J 1 f 

' ^ J/V ax dy ^ /V * 

iT™ 

^ JJ V 2 2 ** ^ Wlm y 2 2 <*y 


(A2) 


(A3) 


From equation (193)> 


r x - Bln g (A4) 

while from equation (624) of reference 25 for the plates considered, in 
general, 


Y„ = 


where 


a 'l* a 2 


= ±A 2 + A(ta/W) 1 / 2 


7lV2 


(a6) 


while A is given hy equation (194). 

For n = 2 .- Considering plate AB in figure 23(a), the conditions 
of simple support at y = 0 require that C-^ = C^ = 0 so that from 
equation (A5) 

Y 2 = C 2 sin a^y + C^ sin a^y (A7) 

At y = b, Yg = 0 so that from equation (A7) 

sinh o£b 

*^4 ~ ~^2 sin a 2 d 


(A8) 
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From table 2 6 of reference 25 the buckling coefficient k for plate AB 
in figure 25 (a) (one edge, simply supported) the other, fixed) is 5 *^ 2 , 
occurring for a half wave length a = {3b = 0.8b so that in equation (A. 6 ) 

tff/N = kn 2 /b 2 = 5.42jr 2 /b 2 


Hence, from equation (A 6 ), with X = it/ 0.8b, one obtains aq = 2.12ic/b 
and otg = l.l65rt/b. Inserting this into equations (A7) and (A 8 ) 


= C 2 sinh (2.12ity/b) + 790 sin (l.l65rty/b) 


Denoting the deflection at y = b/2 as wp^, this becomes 


Y 2 = ^Qsinh ( 2 . 12 :ty/b) + 790 sin (l.l65*y/b)] (A9) 

From equations (A 1 !-) and (A9) one obtains 


[ Y X 2 dy = 0 . 5 bwq m 2 

Jo 

r b 

/ Y 2 2 ay - 0.471^ 

o o 


f Y 1 Y 2 ^ = 0 . 566bwj Lm W2 m 
J Q 


Hence, from equations (A2) and (A3) 9 = 0.752, 7 = 0.778, an <3 

79 = 0.57. 

For n = 00 .- It is sufficiently accurate to calculate 9 and y 
for n = 2 and n = 00 only and to derive from these values the v alue s 
of 9 and 7 for intermediate n values. For infinite values of n 
(fig. 22 (b)) 

Y 1 = w lm = Constant (A10) 

From symmetry (fig. 23(d)) equation (A5) transforms to 


y 2 = ^1 cos k QqY + Cj cos a 2 y 


(All) 
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At y = b/2, Y 2 = 0 so that from, equation (All) 



cos, 


“Cn 


COS 


h ( a l b / 2 ) 

(“2V 2 ) 


(A12) 


From table 26 of reference 25 the buckling coefficient k = 6.97 occurs 
for p = a/b = 0.668, from which, in a way similar to that for n = 2, 
one obtains 


= 2.488rt/b 


ag = 1.503rt/b 


and 

Y 2 ■ ^^jcosh(2.488jty/b) + 54.5 eos(l.303Jty/b)] (A13) 

Hence, from equations (A10) and (A15) 


L 


b/2 


Y X 2 dy - O.^ 2 


,b/2 


f Y 2 2 dy = 0.2065bw 2m 2 
J n 


r b/2 f b/2 

J Q Y 1 Y 2 ^ = w Im J 0 Y 2 ^ = °' 2Tbw M W 2m 


so that from equations (A2) and (A3) 9 = 0.54, 7 = 1.31> and 79 = 0.707* 

For n = 3 »- It is sufficiently accurate to consider the angle 
rotations at B and C (fig. 23(b)) to be zero so that the two outer 
plates AB and CD for n = 3 axe similar to the plates for n = 2 
(fig. 23(a) and the center plate is similar to those for infinite values 
of n (fig. 23(d)). Hence, as an average value 
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For n = 4.- Assuming the angle rotations at B, C, and D to "be zero 
(fig. 23(c)), one obtains, for n = k, 


79 = | (0.57 + 0.707) = 0.639 
Values of 7<p are assembled in table 2. 
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TABLE 1.- COMPARISON OF RESULTS FROM EQUATIONS (46) AND (48) 


h 

cm 

c, 

cm 

n 

b, 

cm 

"■ (W r )/( “=/"<=) 

V 

= 0 

v = 4 x 10”^ cm/kg 

v = 30 x 10 5 cm/kg 

Eq. (46) 

Eq. (48) 

Eq. (46) 

Eq. (48) 

Eq. (46) 

Eq. (48) 

260 

130 

2 

16 

1.11 

1.09 

1.45 

1-37 

1.84 

1.80 

390 

130 

3 

16 

1.27 

1.24 

1.70 

1.63 

2.47 

2.42 

520 

130 

4 

16 

1.47 

1.43 

1.90 

1.85 

2.94 

2.89 

520 

130 

4 

24 

1.24 

1.21 

1.51 

1.46 

2.39 

2.33 

780 

130 

6 

16 

1.92 

I.89 

2.30 

2.26 

3.56 

3.53 

156 

78 

2 

16 

1.12 

1.09 

1.56 

1.48 

1.89 

I.89 

390 

78 

5 

16 

I.69 

1.65 

2.30 

2.26 

3.68 

3.66 

780 

78 

10 

16 

2.94 

2.91 

3-37 

3.34 

5.09 

5.09 


TABLE 2.- VALUES p AND q IN EQUATION (203) 
AND 7<p FROM EQUATIONS (211) AND (221) 


n. 

p 

q 

79 

2 

2.27 

2.45 

0.370 

3 

2.36 

3.58 

.616 

4 

2.37 

4.20 

.639 

00 

2.5 

5 

.707 







2 

2 

3 

3 

4 

1 

4 

4 

CO 

00 

00 

00 

El/(fcN) 

5 

5 

5 

0 

5 

2 

0 

5 

2 

5 

0 

5 

0 

0.4 

0 

0 

0 

0 

0 

0 

0 

0.4 

0.4 

k cr from 

eq.. ( 224 ) . . . 

4.15 

3-73 

4.22 

1.255 

4.25 

2.78 

1.15 

4.25 

2.66 

3.59 

0.71 

k cr from 

ref • • • • • | 

4.20 

3.75 

4.25 

(1.28) 

4.25 

2.79 

1.15 

4.25 

2.66 

3.60 

0.71 
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Figure ] Replacement of "bending and shearing rigidity of sandwich 

column by equivalent spring systems. 



Figure 2.- Long rectangular sandwich plate compressed in X-direction with 
various boundary conditions at loaded and unloaded edges. 
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(t>) Composite case. 

Figure 6.- Clamped sinrply supported sandwich column 
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Figure 8.- Equivalent laterally 
supporting spring system for 
column with batten plates. 



(a) Coupled by 
diagonals . 


\ 


\ 

c 

\ 


\ 


\ 


k 



(b) Coupled by 
diagonals and 
verticals . 


Figure 9.- Latticed columns. 
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Figure 10.- Long sandwich plate compressed in X-direction and clamped 

at unloaded edges. 




Figure 11.- Long sandwich plate with 
isotropic faces and an orthotropic 
core compressed in X-direction 
with arbitrary boundary conditions 
at unloaded edges. 



Figure 14 . - Cross section of 
corrugated sandwich plate. 



Figure l6 . - Buckling of long 

simply supported flange under 
linearly distributed stresses 
that increase with y. 



Figure 15 .- Deflections of Y-strips 
for clamped plates. 



Figure 17-- Buckling of long simply 
supported flange under linearly 
distributed stresses that decrease 
with increasing y. 
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(a) Total deflection. 



(c) Deflection for case 2. 

Figure 20.- Deflection of cross section of I-section. 
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(a) Longitudinal section 
through center of 
buckle . 


(b) Cross section through center 
of buckle. 


CASE / 



(c) Cross section with deflection in case 1. 

CASE Z 



(d) Cross section with deflection in case 2. 
Figure 22.- Buckling deflection of stringer panel. 
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(d) n = 00 . 

Figure 23.- Sketches used in computation of p, q.> and 79 . 
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